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ABSTRACT 

A  separation  algorithm  applicable  to  the  pattern  classification  and 
cluster  analysis  of  n-dimensional  (n  >  2)  data  is  presented.     The  algorithm 
reduces  the  dimensionality  of  the  problem  by  projecting  each  point  into 
a  plane.     This  plane  is  presented  to  the  user  on  a  computer  graphics 
console  screen.    The  operator  picks  a  point  on  the  screen  with  a  lightpen' 
and  chooses  a  "direction  of  movement"  to  achieve  or  increase  separation, 
thereby  causing  an  iteration  of  the  algorithm.     Each  iteration  is  in  fact 
a  reorientation  of  the  plane  into  which  the  data  points  are  projected. 
Iterations  continue  until  satisfactory  separation  is  achieved.    The 
algorithm  is  not  restricted  by  the  dimensionality  of  the  data,  nor  are  any 
distributional  assumptions  required.     Results  from  six  case  studies 
indicate  that  the  algorithm  is  a  useful  tool  for  the  analysis  of  multi- 
dimensional data. 
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I.     INTRODUCTION 

The  separation  algorithm  presented  in  this  paper  is  a  development 
of  an  idea  of  Professor  Thomas  M.  Cover  of  Stanford  University.    It  pro- 
vides a  tool  for  classifying  multidimensional  data  sets.    Two  classification 
methods  can  be  handled  by  the  program  which  implements  the  algorithm; 
these  two  methods  are  pattern  classification  and  cluster  analysis. 
Meisel  [4]  discusses  and  contrasts  these  two  processes.     He  describes 
pattern  classification  as  the  "process  of  developing  on  the  basis  of  a 
finite  set  of  labelled  samples,  a  decision  rule  with  which  we  can  classify 
a  point  in  the  pattern  space  corresponding  to  an  unlabelled  sample.  " 
Cluster  analysis  on  the  other  hand  is  described  as  a  process  for  "describing 
and  locating  the  discernible  subsets  of  which  a  set  of  sample  points  of 
reasonable  complexity  is  usually  formed." 

Meisel  notes  that  the  difference  between  pattern  classification  and 
cluster  analysis  is  that  the  samples  for  the  former  are  labelled,  whereas 
in  the  latter  they  are  unlabelled.    Thus  it  appears  that  the  question  of 
whether  a  particular  problem  should  be  approached  as  a  pattern  classifi- 
cation or  cluster  analysis  problem  depends  on  the  use  to  be  made  of  the 
results  and  on  the  nature  of  the  data  available.    In  the  pattern  classifi- 
cation problem  the  initial  data  set  is  divided  into  two  (or  more)  subsets 
which  are  identified.    Each  member  of  the  data  set  is  labelled  as  being 
included  in  exactly  one  of  these  subsets.    The  decision  rule  developed 
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typically  takes  the  form  of  a  function  which  has  an  argument  in  the  form 
of  a  data  point.    When  input  with  the  values  of  elements  of  different 
subsets,  the  function  assumes  values  in  disjoint  ranges.    For  example  if 
the  elements  of  a  two  subset  problem  are  identified  as  X  and  Y  respectively, 
and  the  discriminant  function  is  f ,  possibly  f(X)  would  be  negative  for  all  X; 
and  f(Y)  would  be  positive  for  all  Y.     The  discriminant  function  is  then 
used  to  classify  unlabelled  points  Z  as  being  in  either  the  X  or  the  Y  sub- 
set. 

In  a  cluster  analysis  problem,  the  initial  data  set  has  no  previously 
defined  subdivisions.     The  purpose  of  the  analysis  is  to  decompose  the 
given  set  into  a  group  of  subsets.    The  number  of  subsets  to  be  found 
can  be  specified  in  advance  or  unknown.     Utilizing  similarities  among  the 
elements  in  the  original  data  set,  assignment  to  the  subsets  is  made. 

Both  pattern  classification  and  cluster  analysis  techniques  can  be 
described  as  being  direct  or  indirect.    In  an  indirect  process,  some 
criterion  function  is  used  in  the  construction  of  the  discriminant  function 
or  is  used  to  define  the  quality  of  clustering  desired.    A  direct  process 
has  no  such  criterion  function;  the  formulation  of  the  discriminant  function 
(or  clusters)  is  accomplished  by  making  immediate  use  of  the  data  points. 
The  algorithm  developed  in  this  paper  is  direct  in  nature. 

It  may  be  that  for  a  particular  problem,  both  pattern  classification 
and  cluster  analysis  techniques  are  appropriate.    For  example  in  a  pattern 
classification  problem  with  two  subsets  forming  the  data  set,  the  subsets 
might  be  inspected  separately  using  cluster  analysis.    Alternatively 


after  successful  usage  of  clustering  techniques,  a  pattern  classification 
algorithm  might  be  employed  to  define  discriminant  functions  for  the 
clusters . 

The  separation  algorithm  presented  in  this  thesis  is  implemented 
by  a  computer  program  which  utilizes  the  graphics  edit  mode  of  operation 
provided  by  the  XDS93  00  digital  computer  combined  with  the  Adage  AGT-10 
graphics  console.    The  man-machine  interface  is  emphasized.     The  basic 
concept  of  the  algorithm  becomes  apparent  when  a  comparison  is  made 
with  the  concept  of  Fisher's  Linear  Discriminant.    Duda  and  Hart  [3] 
present  Fisher's  Linear  Discriminant  as  a  method  for  reducing  the  dimen- 
sionality of  a  pattern  recognition  problem.     By  projecting  the  points  of 
the  n-dimensional  pattern  space  onto  a  straight  line  in  the  pattern  space, 
the  dimensionality  of  the  problem  is  reduced  from  n  to  one.     They  further 
point  out  that  the  projections  of  the  points  may  fall  on  the  line  in  a  very 
confused  manner.     By  changing  the  position  of  the  line  in  the  n-space, 
an  orientation  is  sought  in  which  the  projections  of  the  data  points  will 
display  separation.    The  achievement  of  this  orientation,  if  possible,  is 
the  goal  of  the  discriminant  analysis. 

By  way  of  comparision,  the  algorithm  presented  here  also  represents 
a  dimension  reduction  approach;  however,  the  projection-of  the  data 
points  is  made  onto  a  plane  vice  a  line.     The  dimension  reduction  is  then 
from  n  to  two.    If  n  is  less  than  or  equal  to  two,  this  algorithm  is  not 
appropriate;  more  straightforward  methods  employing  scatter  plots,  etc.  , 
can  be  employed.    As  in  the  Fisher  method,  the  algorithm  proceeds  as  a 


series  of  reorientations  of  this  plane  in  the  n-space.     Following  each  new 
orientation  the  program  presents  the  projections  on  the  screen  of  the 
graphics  console.     In  a  pattern  recognition  problem  with  two  subsets, 
each  point's  projection  is  represented  by  a  video  "Y"  or  "X".     In  a  cluster 
analysis  each  projection  is  presented  as  a  "Y"  .    With  the  XDS9300/AGT-10 
combination  in  the  graphics  edit  mode,  the  user  decides  whether  or  not 
to  make  an  iteration,  and  thus  reorient  the  plane.    The  decisions  if  and 
how  to  iterate  are  based  on  the  situation  presented  on  the  console  screen. 
The  importance  of  the  man-machine  interface  should  be  clear.    The  iterative 
process  will  be  described  in  greater  detail  later.    Note  that  the  decision 
of  how  to  iterate  and  when  to  terminate  rests  with  the  user.    He  determines 
when  adequate  separation  has  been  attained. 


II.    A  PREVIOUS  COMPUTER-GRAPHICS  APPLICATION 

This  discussion  of  a  report  by  Herman  Chemoff  entitled  "THE  USE 
OF  FACES  TO  REPRESENT  POINTS  IN  n-DIMENSIONAL  SPACE  GRAPHICALLY" 

1    is  presented  as  an  example  of  the  type  of  man-machine  graphics 
approach  to  multidimensional  data  analysis  previously  proposed.     Chernoff 
uses  each  n-dimensional  vector  corresponding  to  a  data  point  to  mathe- 
matically determine  the  shape  and  features  of  a  human  face  which  is  output 
by  a  CALCOMP  plotter.     Chernoff  s  hypothesis  is  that  points  belonging  to 
the  same  class  or  cluster  will  be  identified  by  the  user  due  to  the  similarity 
of  the  faces  they  generate.    The  user  constructs  the  clusters  by  grouping 
the  faces  into  sets  which  show  common  features.    As  an  illustrative 
example,  Chernoff  applies  his  method  to  a  data  set  of  eighty  eight- 
dimensional  measurements  of  nummulitid  specimens  from  the  Eocene 
Yellow  Limestone  Formation  of  Northwestern  Jamaica.     (This  data  set  will 
also  be  used  later  as  an  illustrative  example  for  this  paper.)    Chernoff's 
results  indicated  the  existence  of  three  clusters  in  the  data.    These  results 
correspond  exactly  to  the  results  obtained  by  Wright  and  Switzer  [6]  who 
employed  a  different  cluster  analysis  scheme  to  the  same  data  set. 

In  constructing  faces  Chernoff  can  handle  up  to  eighteen  dimen- 
sional data  to  determine  such  features  as  the  shape  and  spacing  of  the 
eyes,  size  and  curvature  of  the  mouth,  and  the  shape  of  the  head. 
Clustering  obtained  by  Chernoff's  method  may  depend  on  which  components 
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of  the  data  vectors  control  the  various  facial  features;  it  is  not  clear 
how  this  assignment  should  be  made  before  the  problem  is  started. 
Chemoff  dismisses  any  criticism  of  the  eighteen  dimension  limitation, 
saying  that  any  further  increase  in  dimensionality  could  be  overcome  by 
adding  features  to  the  faces.     This  is  an  interesting  point;  it  would  seem 
that  if  too  many  features  were  added  to  the  faces,  the  decision  process 
would  be  more  complicated.     Chernoff  notes  that  "when  the  eyes  are  very 
small,  the  pupils  become  hard  to  detect."    Also  he  notes  that  some  infor- 
mation is  lost  when  the  two  ellipses  which  form  the  head  are  not 
differentiable  and  the  face  is  circular.    Certainly,  within  any  moderate 
range  of  dimensional  values,  and  as  noted  by  Chernoff,  this  method 
"provides  a  promising  approach  for  a  first  look  at  multivariate  data  which 
is  effective  in  revealing  rather  complex  relations  not  always  visible  from 
simple  correlations  based  on  two-dimensional  linear  theories.  " 


III.     DISCUSSION  OF  THE  ALGORITHM 

As  indicated  in  Section  I  the  separation  algorithm  developed  in  this 
paper  is  applicable  to  both  the  pattern  classification  and  the  cluster 
analysis  problems.    The  concept  of  operation  of  the  algorithm  in  each  of 
these  applications  is  the  subject  of  the  following  paragraphs. 

For  the  pattern  classification  problem  the  analysis  is  conducted 
on  an  n-dimensional  data  set  which  is  divided  into  two  subsets.    These 
subsets  might  represent  two  classes  of  objects,  P     and  P    .    Each  object 
would  be  characterized  by  an  n-vector  (or  pattern)  of  measurements.    Each 
element  of  the  pattern  is  referred  to  in  the  literature  as  a  feature.     The 

goal  of  the  pattern  classification  analysis  is  to  determine  a  function 

_        i  _ 

(called  the  discriminant  function)  f:R— >R    ,  where  R    is  real  n-space, 

which  can  be  used  to  assign  an  unlabelled  pattern  Z  to  either  P    or  P    . 

A  scheme  for  assigning  Z  might  be: 

f(Z)  <0=>  Z«P 

f(Z)>0=^Z€P2 
If  a  discriminant  function  of  reasonably  low  complexity  can  be  determined 
for  the  two  classes,  the  classes  are  separable;  if  no  such  function  is 
found,  the  classes  are  inseparable.    Since  the  discriminant  function  is 
constructed  using  samples  from  each  of  the  classes,  it  must  be  determined 
that  the  separation  between  the  classes  is  sufficient  to  warrant  the  use 
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of  the  discriminant  function  in  making  class  assignments.     For  example, 

for  two  classes  to  be  linearly  separable,  the  convex  hulls  of  the  two  sets 

of  projections  should  be  disjoint. 

For  purposes  of  illustration  and  discussion,  suppose  the  initial 

data  set  has  N  members  in  P     and  M  members  in  P    .    The  P1   subset  is 

read  by  the  computer  as  Y  vectors;  thus  the  program  reads  N  vectors  Y. , 

i  =  1  ,N  .     Similarly  the  P9  subset  is  read  as  a  set  of  X  vectors  ,  X. , 

j  =  1,M.     The  dimension  n  of  these  vectors  should  be  greater  than  two. 

(If  n  is  two  or  less  the  problem  can  be  handled  graphically.)    The  algorithm 

proceeds  to  the  dimension  reduction  task  as  follows.    The  user  chooses  two 

non-zero  vectors  A  and  B  in  Euclidean  n-space.    A  linear  transformation 

is  established  as  TR    „(Z)  =  (A*  Z,B- Z)  where  A-Z  =      z_   a. z.  is  the  inner 

A,B  t=l      1   1 

product  and  Z£R    .    This  transformation  which  is  a  projection  of  the  point 
Z  into  the  plane  defined  by  the  vectors  A  and  B  is  the  algorithm's  vehicle 
for  dimension  reduction.    The  projection  is  not  generally  orthogonal.     The 
transformation  is  applied  to  each  member  of  P    and  P  ;  the  two  sets  of 
coordinate  pairs  are  then  scaled.    In  order  to  be  presented  on  the  screen 
of  an  AG T  console,  a  coordinate  pair  (x,y)  must  have  values  -l^x,yfi  1; 
it  is  to  obtain  pairs  within  this  range  that  the  scaling  is  done.     The 
scaling  is  accomplished  by  searching  through  both  sets  of  projection 
pairs  to  find  the  maximum  absolute  value  of  the  first  coordinates  (called 
TMAX)  and  the  maximum  absolute  value  of  the  second  coordinates  (called 
SMAX).     Then  the  typical  (A-X,B-X)  becomes  (A  'X/TMAX,B  •  X/SMAX) , 
and  the  typical  (A-Y,B-Y)  becomes  (A*  Y/TMAX,B  •  Y/SMAX) .     Thus  the 
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range  of  values  for  all  pairs  of  the  problem  is  as  required  by  the  AGT 
console.    The  resulting  A-B  plane  is  graphically  presented  to  the  user  on 
the  AGT  screen.    It  is  the  user's  decision  based  on  information  gained 
from  the  screen  to  decide  how  to  proceed.     In  all  but  the  most  extraordinary 
case  the  initial  picture  will  show  some  intermingling  of  the  X's  and  Y's. 
The  heuristic  argument  which  follows  provides  a  reasonable  guide  for 
proceeding. 

The  user  is  seeking  a  presentation  which  depicts  the  X's  and  Y's 
in  disjoint  concentrations.     Frequently  the  initial  view  of  the  A-B  plane 
will  yield  an  indication  of  where  these  concentrations  will  develop.    If 
this  is  not  the  case,  and  provided  there  is  a  solution,    such  an  indication 
should  appear  after  a  sequence  of  "random"  iterations.     The  user  then 
picks  the  X  or  Y  which  is  deepest  in  the  opposing  concentration  and 
"moves"  it  toward  its  proper  concentration.    The  user  implements  these 
decisions  with  the  box  of  sixteen  function  switches  and  the  lightpen 
attached  to  the  console.    The  positions  of  the  function  switches  and  the 
lightpen  relative  to  the  graphic  console  are  indicated  in  Figure  III—  1 . 
The  user  makes  his  choice  by  depressing  the  function  switch  marked 
"PICK  X"  or  "PICK  Y"  and  placing  the  tip  of  the  lightpen  over  the  X  or  Y. 
Figure  III-2  shows  a  user  at  the  console  ready  to  pick  a  point  with  the 
lightpen.    The  point  seen  by  the  lightpen  partially  vanishes;  that  is,  a 
Y  might  become  a  v  and  an  X,  a  Y  .    When  part  of  the  desired  point 
disappears,  the  user  completes  the  choice  by  depressing  the  switches 
marked  "COMPLETE  CHOICE"  and  "END  EDIT".    Ready  to  move,  the 
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Figure  III- 1 .    AGT  CONSOLE 


Figure  III-2  .     OPERATOR  SITTING  AT  AGT  CONSOLE 
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user  depresses  one  of  four  switches  (LEFT,  RIGHT,  UP,  or  DOWN)  to 
indicate  the  desired  direction  of  motion.     Figures  III— 3  and  III-4  give  an 
indication  of  what  would  occur  on  two  iterations.     In  Figure  III—  3  the  user 
picks  the  Y  indicated  and  iterates  with  a  move  to  the  right.     Then  in 
Figure  III-4  he  picks  an  X  and  iterates  with  a  move  to  the  left.     The 
iteration  step  consists  of  the  modification  of  the  A  and  B  vectors,  the  re- 
computation  and  rescaling  of  the  projections,  and  finally  the  representation 
of  these  projections  on  the  AGT  screen.    The  modified  values  of  the  A 
and  B  vectors  depend  on  the  point  picked  and  the  direction  of  motion  chosen 
If  the  point  picked  is  the  projection  of  data  point  Z,  the  new  A  and  B 
vectors  will  be  as  follows: 

Direction  of  Move  NEW  A  NEW  B 

LEFT  A-Z  B 

RIGHT  A+Z  B 

UP  A  B+Z 

DOWN  A  B-Z 

The  new  A  and  B  vectors  are  used  in  the  recomputation  of  the  T  transfor- 
mations . 

Having  achieved  the  desired  separation,  the  user  causes  the  line 
printer  (XDS9300  output  device)  to  print  a  six  by  six  inch  reproduction 
of  the  projections  in  the  A-B  plane.     Using  this  printout  he  may  fit  a 
suitable  discriminant  function  between  the  two  concentrations.    In 
Figure  III- 5  an  output  with  the  discriminant  function  is  represented.     The 
discriminant  function  f(Z)  in  this  case  is  a  straight  line,  and  its  equation 
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Figure  III- 3.  Iteration  by  picking  a  Y  and  moving  it  right 
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Figure  III-4.    Iteration  by  picking  an  X  and  moving  it  left 
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Function:  f=0 
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Figure  III- 5.  Output  with  discriminant  function 
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would  be: 

f(Z)  =    OC  (A-  Z/TMAX)  +     p  (B  •  Z/SMAX)  +  tf 
where     PC   ,     P  ,  and    ft  are  scalars.    The  vectors  A  and  B  and  the  scaling 
factors  TMAX  and  SMAX  have  values  corresponding  to  the  printout.    With 
the  specification  of  the  discriminant  function,  the  pattern  classification 
is  completed. 

The  employment  of  the  separation  algorithm  for  cluster  analysis  is 
very  similar  to  that  described  for  pattern  classification.     For  a  cluster 
analysis  problem  the  initial  set  of  data  points  is  not  broken  down  into 
subclasses;  the  program  reads  all  patterns  as  Y  vectors.    The  iterative 
operation  of  the  cluster  separation  is  the  same  as  the  pattern  classifi- 
cation except  that  there  are  no  X's  present.    When  the  user  considers 
that  he  has  displayed  a  satisfactory  set  of  clusters,  he  can  have  the  line 
printer  output  the  six  by  six  inch  picture.    The  output  also  provides  a 
listing  by  point  of  the  projection  coordinates.     Thus  the  cluster  member- 
ship of  each  point  can  be  determined.    As  in  the  pattern  classification 
problem,  no  criterion  function  is  provided  to  direct  the  procedure  or  to 
signal  termination.    It  rests  with  the  user  to  make  value  judgments  to 
conduct  the  analysis.    There  is  no  guarantee  of  a  unique  clustering 
scheme  for  a  given  data  set;  comparisons  between  different  schemes 
should  provide  a  great  deal  of  information. 

As  a  final  comment  on  the  verbal  description  of  the  two  types  of 
operations  just  presented,  it  should  be  noted  that  in  reality  the  program 
does  not  "move"  any  points.    Instead,  it  is  reorienting  the  plane  into 
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which  the  points  are  projected  by  modifying  the  A  and  B  vectors.     Separ- 
ation is  achieved  when  the  A-B  plane  is  properly  oriented. 

Mathematically  the  progress  of  the  algorithm  can  be  traced  through 
a  series  of  steps  as  shown  below.    Suppose  the  initial  data  set  consists 
of  k-dimensional  vectors  Y.  ,  i=l,N  andX.,  j=l,M.     In  the  cluster  analysis 
problem  M=0.    With  the  non-zero  k- vectors  A  and  B  specified,  the 
algorithm  proceeds: 

1)  Center  the  data  sets  by  calculating  a  k-vector 

C=(1/M+N)(Z_    Y.+  LX,L  and  use  C  to  center  the  X's 
i*i      i    j-t     J 

and  the  Y's  as  follows: 
X'=X.-C,  j=l,M 

y:=y:-c,  i=i,N 

11 

2)  Calculate  the  projection  T         for  all  X  and  Y: 

A,  B 

T.    _(X!)  =  (A'X'.,B-X:)  j=l,M 

T.  r(y:)  =  (a-y:,b-y:)  i=i,n 

J\ ,  D        1  1  1 

3)  Calculate  the  scaling  parameters  TMAX  and  SMAX 

TMAX  =  max(max|A-X'.  |  fmax[A-Y:  |) 


i=l,N 
J=1,M 


SMAX  =  max(max|A'X'.  |  ,  max[B-Y:|   ) 


4)         Scale  the  T         projections 


TA    _(X:)  =  (A'X'/TMAX,  B'X'./SMAX)  j=l,M 

«,D        J  J  J 

T.    D(Y!)  =  (A'Y:/TMAX,  B-Y!/SMAX)  i=l,N 
A,B     i                 i  i 
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5)         At  this  point  M  video  X's  and  N  Y's  appear  on  the  AGT  screen, 
if  a  Y  is  picked  and  a  right  move  is  selected,  a  new  A=A+Y 
is  calculated.     Thus  at  this  point,  either  a  new  A  vector  or 
a  new  B  vector  is  found  and  control  shifts  to  step  2  . 

The  five  steps  listed  above  describe  the  activity  in  the  program  at  each 

iteration. 
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IV.     MOTIVATION  OF  THE  ALGORITHM 

The  attraction  of  the  approach  of  this  algorithm  for  separating  data 
into  its  classes  has  a  dual  basis.    Study  of  the  iterative  operation  of  the 
algorithm  reveals  as  the  first  part  of  this  basis  the  manner  in  which  the 
correlation  within  data  classes  is  exploited  in  seeking  separation.     The 
second  grounds  for  confidence  is  the  similarity  of  the  algorithm's  method 
to  that  proposed  in  the  Perceptron  Convergence  Theorem.    This  similarity 
will  be  utilized  to  present  a  convergence  proof  for  the  algorithm .    Each 
of  these  supporting  concepts  will  be  discussed  in  the  following  paragraphs 

If  the  data  set  to  be  analyzed  via  the  algorithm  is  separable,  the 
set  can  be  considered  to  be  made  up  of  a  number  of  subsets.     For  illus- 
trative purposes,  suppose  the  data  set  consists  of  two  subsets.    The 
assumption  is  made  that  there  is  a  higher  correlative  measure  within 
these  subsets  than  between  them;  it  is  this  difference  which  has  value 
in  the  attainment  of  separation.    If  the  two  subsets  are  labelled  X  and  Y, 
then  without  considering  the  scaling  operation,  it  is  clear  from  the  dis- 
cussion of  the  last  section  that  the  algorithm  maps  these  sets  onto  a 
plane  as 

(A'X.,B-X.)    forX.e   X 
J  J  J 

(A-Y.  ,B-Y.)    forY.  €    Y 
11  1 

Now  suppose  an  iteration  is  made  by  picking  an  X  and  moving  to  the  right. 
It  has  been  seen  that  the  new  A  is  in  fact  A+X,  .    The  new  mappings  are 
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((A+X,  )'X.,B-X.)  =  (A-X  +X_  -  X.,B'X.),  for  X.  €    X 
k      ]  j  j      k      j  j  j 

((A+X,)«Y  ,B'Y.)  =  (A-Y.+X^  Y.  ,B'Y.),  for  Y    €    Y 

The  difference  in  the  mappings  (i.e.  ,  the  separation  gained)  is  seen  to 

be  in  the  first  coordinates  as  the  difference  between  X,  •  X.  and  X  •  Y. . 

k      j  k      1 

Standard  correlation  concepts  indicate  that  if  the  X  set  is  highly  corre- 
lated, while  the  correlation  between  the  X  and  Y  sets  is  low,  then  the 
product  X*  X.  will  be  greater  than  X»  Y.  .    Clearly  if  some  Y     is  picked 
the  difference  in  the  mappings  caused  by  the  iteration  will  be  the  difference 

between  Y  *  Y,  and  Y  *  X..    In  general  the  degree  of  separation  gained 
m      1  m      j 

on  a  single  iteration  will  tend  to  be  small;  however,  it  is  the  goal  of  the 
algorithm  to  combine  the  incremental  gains  in  separation  of  a  sequence 
of  iterations  to  obtain  the  desired  distinction  between  classes. 

To  perceive  the  similarity  of  the  algorithm  to  the  method  of  the 
Perceptron  Convergence  Theorem,  one  has  only  to  consider  the  following 
version  of  this  work  as  presented  by  Minsky  and  Papert  [5]  : 
"The  Perceptron  Convergence  Theorem" 

Consider  the  following  program  which  the  vector 
notation  A-<£  is  2_°*^'    (i-n  place  of  our  usual  notation) 
START  Choose  any  value  for  A 
TEST  Choose  an  X  from  F  U  F~ 

If  X  £  F+  and  A  -  i  >  0    go  to  TEST 
If  X€.  F+  and  A-#  £0    go  to  ADD 
If  X  €  F~  and  A  •  $  <■  0    go  to  TEST 
If  X  £  F"  and  A'$£0    go  to  SUBTRACT 
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ADD  Replace  A  by  A  +  #  (X) 

Go  to  TEST 
SUBTRACT  Replace  A  by  A  -  $  (X) 
Go  to  TEST 
We  assume  until  further  notice  that  there  exists  a  vector  A*  with 
the  property  that  if  X  €  F+  then  A*-  *  (X)  >  0  and  if  X  £  F~  then 
A**  $    (X)  <  0.     The  perceptron  convergence  theorem  then  states 
that  whatever  choice  function  is  used  in  TEST,  the  vector  will 
be  changed  only  a  finite  number  of  times .  " 
The  variable  A  of  the  above  theorem  corresponds  to  the  A  of  the  present 
algorithm;  there  is  no  B  vector  in  the  Perceptron  theorem  approach.    The 
Perceptron  method  can  be  considered  to  be  a  more  direct  method  than 
that  presented  for  the  algorithm  here  since  it  specifies  that  only  mis- 
classified  points  will  be  operated  on.     Naturally  this  approach  is 
allowable  but  not  required  in  the  present  work.    If  it  is  assumed  that, 
in  employing  the  algorithm,  the  user  conforms  to  the  modus  operandi  of 
the  Perceptron  Convergence  Theorem,  it  is  possible  to  prove  that  if  the 
data  is  linearly  separable  the  algorithm  will  attain  separation  within  a 
finite  number  of  steps.     The  following  proof  is  a  modification  of  the 
method  employed  by  Duda  and  Hart  [3]  ;  it  shows  that,  provided  there 
exists  an  orientation  of  the  A-B  plane  in  which  the  X  subset  is  linearly 
separable  from  the  Y  subset,  this  orientation  can  be  attained  with  a 
finite  number  of  iterations.    Suppose  the  situation  in  the  A-B  plane  is 
as  indicated  in  Figure  IV- 1.    Let  b=B/SMAX  and  a=A/TMAX.    Then 
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Origin  v.                ^\ 

\  Y'S    ) 

\  X'S    ) 

"" 

Figure  IV- 1.    Orientation  of  the  A-B  plane  in  which  the 
X  and  Y  sets  are  linearly  separated 
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b»Z  -  m(a*  Z)  +  w    =  f(Z) 
o 

where  m  is  the  slope  of  f,  w     is  the  y  intercept,  and  Z  is  an  X  or  Y.    Then 

f(Z)  =  (b-ma)'Z  +  w 

o 

To  simplify  this  expression  the  following  n+1  dimensional  (where  n  is  the 

dimension  of  A,B,  and  Z)  vectors  are  defined: 

m      ,b-ma>  .    _       ,Z. 

T  =  (          )  and   S  =  (J 

w  i 
o 

Then  clearly  we  can  write 

f'(S)  =  T'S 
Then  f'(S)  is  a  hyperplane  which  passes  through  the  origin  of  the  n+1 

dimensional  "S"  space.    As  Duda  and  Hart  point  out,  the  addition  of  the 

st 
n+1      dimension  as  a  one  to  the  n-dimensional  Z  preserves  all  distance 

relationships  among  the  points  of  the  original  data  set.    The  S  vectors 

are  all  in  the  same  n-dimensional  subspace  of  the  n+1  space.    This 

translation  into  n+1  space  reduces  the  problem  from  examining  the  vectors 

a  and  b  and  the  scalars  w    and  m  to  the  problem  of  inspecting  the  vector 

X  Y 

T.    From  Figure  IV- 1  it  is  apparent  that  f  (   )<0  for  all  X,  and  f  (  )>0  for 

all  Y. 

Now  defining  a  variable  U  such  that: 
U  =  f  *)  for  all  X 

U  =  m     for  all  Y 

the  problem  is  confined  to  discussing  f'(U)  =  T»U  >  0  for  all  U.     The 

Y 
behavior  of  the  user  is  such  that  if  he  picks  a  Y  such  that  f  (  )  <  0,  he 

causes  the  program  to  iterate  with  T  =  T+Y.    Similarly  if  he  picks  X  such 
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that  f  (,)  >  0,  he  iterates  with  T  =  T-X.     In  either  case  the  iteration  can 

be  considered  as  picking  a  U  and  iterating  with  T  =  T+U.    The  iterations 

are  conducted  only  when  the  function  as  formed  misclassifies  the  U  picked 

Let  T.    be  the  value  of  the  T  vector  on  the  k      iteration;  U.    is  the 
k  k 

th  A 

point  picked  for  the  k      iteration;  and  e  is  a  positive  scalar.    Also  T  is 

the  "solution"  vector  T.    Then 

(Tk+1-eT)  =  (Tk-eT)+Uk  (a) 


:k+1-eT||2=    ||Tk-eT||2+2(Tk-eT>Uk+     Ifu^ll   2 


However  T.  •  U.  <  0 .     Thus 
k       k 

l|Tk+1-eT||  2  £   ||Tk-ef  ||  2  -2e$-Uk   +  ||U]cll  2  (c) 

A 

Since  T«U,    is  positive,  the  second  term  on  the  right  of  equation  (c)  will 

dominate  the  third  if  e  is  sufficiently  large.    In  particular  let  r=max  U. 

i 

A 

and  s=min  T-U.    for  i=l ,  .  .  .  ,m  +m?  (assume  there  are  m  X's ,  and  m9Y's) 

l 
Then 

Cn  2   s       ii  ^,,2 


||T        -elll  "i  .||T-ejJ|  2-2es  +  /  (d) 


Now  if  e=r  /s 

Tk+1-eT||2  ^     ||Vef||2-r2  (e) 

A  2 

Thus  the  squared  distance  from  T,    to  eT  is  reduced  by  at  least  r    at 

each  iteration.    After  K  iterations 

||Tk+1-eT||2  ^     ||TreT||2-Kr2  (f) 

Since  the  squared  distance  can  not  become  negative,  it  follows  that 

the  sequence  must  terminate  after  no  more  than  K     iterations ,  where 

o 

Kq=   ||Tl-eTl|  2/r2 
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Thus  the  number  K    provides  an  upper  bound  on  the  number  of 
iterations  required  to  transform  the  initial  vector  T     to  the  solution 

A 

vector  T. 
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V.     DISCUSSION  OF  THE  COMPUTER-GRAPHICS  PROGRAM 

The  construction  of  the  program  formulated  to  implement  the  separa- 
tion algorithm  will  be  discussed  in  this  section,  primarily  to  demonstrate 
what  operations  are  performed  on  the  data  during  the  execution  of  the 
program .    Instructions  for  the  use  of  the  program  are  provided  in  Appendix 
II. 

The  program  is  written  in  FORTRAN  IV  as  modified  for  the  XDS9300 
computer  installation  at  the  Naval  Postgraduate  School.    It  employs  the 
library  graphics  program  "MAD"  to  allow  use  of  the  AGT  graphics  console 
with  function  switches  and  lightpen  to  alter  the  graphics  display  and  to 
provide  printed  output  via  the  line  printer  of  the  display.    The  main 
program  can  be  considered  to  consist  of  four  parts: 

1)  Set  Up 

2)  Projection  Calculation  And  Display 

3)  Display  Modification 

4)  Output 

In  the  set  up  portion  storage  is  established  in  the  9300  to  handle 
arrays  based  on  the  dimension  and  size  of  the  set  or  sets,  if  there  are 
two,  of  data.    The  data  is  read  into  the  computer  via  two  formatted  READ 
statements.    The  first  data  set  is  labelled  Y  and  the  second  is  X  (for  the 
pattern  classification  problem).    If  a  cluster  analysis  problem  is  to  be 
conducted  the  variable  M  (number  of  members  in  the  second  set)  is  set 
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to  zero.    The  program  then  reads  the  data  in  as  one  group  of  Y's.    The 
only  modification  made  directly  to  the  data  follows.    The  data  as  a  whole 
is  centered.     For  example  in  three  dimensions,  three  averages  DAV(l), 
DAV(2),  and  DAV(3)  are  computed  over  the  elements  of  the  members  of 
the  population  of  the  X's  and  Y's.    Then  each  Y=(y    ,y„  ,y,J  is  modified 
to  (y   -DAV(l),y2-DAV(2)/y3-DAV(3))  and  similarly  X=(xx  ,x2  ,x3)  is 
translated  to  (x  -DAV(l),x2-DAV(2),x3-DAV(3)). 

The  final  inputs  of  the  set  up  section  are  the  initial  values  of  the 
A  and  B  vectors .    The  program  then  proceeds  into  the  calculation  of  the 
projection  coordinate  pairs  of  each  Y  and  X;  the  second  part  of  the  program 
now  starts . 

Through  a  succession  of  DO-loops  the  projections  of  the  Y's  and 
X's  are  calculated.    The  AGT  screen  is  set  up  to  display  two  dimensional 
points  (G,H),  where  -1—  G,H  —  1.    This  necessitates  that  the  projections 
be  appropriately  scaled  at  each  iteration.    The  scaling  is  accomplished 
by  searching  through  all  pairs  for  the  maximum  absolute  values  of  each 
of  the  two  coordinates.    The  first  is  called  TMAX,  the  second,  SMAX. 
Then  if  a  point  Y  maps  into  a  projected  point  (YY1,YY2),  the  point  dis- 
played on  the  AGT  screen  will  have  coordinates  (YY1/TMAX,YY2/SMAX) . 
The  program  then  proceeds  to  display  the  set  of  the  projections  generated 
from  the  Y's  followed  by  the  projections  generated  from  the  X's.    With 
the  presentation  of  these  points,  the  user  using  the  function  switches 
and  the  lightpen  can  pick  an  X  or  Y  and  cause  an  iteration  by  picking  the 
direction  of  shift  (LEFT,  RIGHT,  UP,  or  DOWN);  this  occurs  in  the  third 

portion  of  the  program  . 
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If  the  operator  chooses  to  make  another  iteration  by  picking  an  X 
or  Y  on  the  screen  and  specifying  a  direction,  the  program  identifies  the 
data  vector  corresponding  to  the  chosen  point  and  using  the  centralized 
version  of  this  data  point,  calculates  the  new  A  or  B  vector.    A  loop  is 
then  made  back  to  the  beginning  of  the  second  section.     There  are  options 
in  addition  to  making  an  iteration  open  to  the  user.    He  can  cause  the 
program  to  loop  back  to  the  first  section  for  inputting  new  values  for  the 
A  and/or  B  vectors.     He  does  this  by  depressing  the  switch  "NEW  A&B". 
By  pushing  the  function  switch  "REGRESS",  he  can  return  to  the  display 
previous  to  the  present.     Finally  by  pushing  "OUTPUT"  he  can  cause  the 
program  to  proceed  to  section  four. 

In  section  four  the  program  fills  a  sixty-one  by  thirty- seven 
dimensional  array  with  symbols  corresponding  to  the  present  AGT  view. 
This  array  prints  out  as  a  six  by  six  inch  square.    The  square  is  outlined 
by  asterisks  except  for  four  zeros  indicating  the  intercepts  of  the  coordinate 
axes  with  the  edge  of  the  square.    The  points  within  the  square  appear  as 
an  X,  0,  J,  8,  or  E.     The  X  and  0  correspond  respectively  to  the  position 
of  an  X  or  Y  projection  in  the  array  cell.    An  8  corresponds  to  two  or  more 
Y's  in  a  cell;  an  E  implies  two  or  more  X's.    A  J  indicates  that  a  combina- 
tion of  X's  and  Y's  occupies  the  cell  in  question.    When  the  printout  is 
complete,  control  is  shifted  to  section  three  to  allow  the  user  to  decide 
what  to  do  next. 

It  is  clear  that  the  only  modification  to  the  data  is  the  centralization 
to  deviations  about  the  means  performed  in  section  one.    The  projection 


29 


of  an  n-dimensional  point,  for  example,  Y,  has  coordinates 
(A-(Y-Y)/TMAX,B*  (Y-Y)/SMAX) .      As  mentioned  earlier,  this  section  of 
the  thesis  is  only  intended  to  describe  the  operation  of  the  program  and 
to  depict  the  relation  between  an  initial  data  point  and  its  eventual 
projection. 

The  compilation  time  for  this  program  is  approximately  eight  minutes 
Of  course,  the  total  running  time  for  the  program  varies  with  the  dimen- 
sionality of  the  problem,  the  size  of  the  data  set(s),  and  with  the  number 
of  iterations  required  by  the  user.    The  running  time  per  iteration  in  a 
six  dimensional  problem  with  two  hundred  and  eighty  data  points  was 
found  to  be  approximately  twenty  seconds.    A  user  can  expect  to  spend 
between  thirty  and  forty-five  minutes  making  iterations  on  a  problem  of 
this  magnitude.    If  there  is  no  reasonably  clear-cut  separation  between 
the  sets  of  a  pattern  classification  problem  or  between  the  clusters  of  a 
cluster  analysis,  the  running  time  may  be  considerably  longer. 


30 


VI.     EXAMPLES  AND  APPLICATIONS 

As  a  means  of  verifying  the  proper  operation  of  the  separation 
algorithm  in  the  graphics  program,  four  data  sets  were  generated  for 
pattern  classification  analysis  as  follows: 

1)  Two  three-dimensional  cubes 

2)  Two  six-dimensional  cubes 

3)  Four  six-dimensional  cubes 

4)  Four  dimensional  sphere  within  a  four-dimensional  shell 
Each  of  these  data  sets  possessed  geometric  qualities  which  could 

easily  be  perceived,  thus  allowing  concrete  evaluation  of  the  operation 
of  the  program  on  multidimensional  data.    Before  addressing  each  data 
set  individually,  a  few  general  comments  are  appropriate.    Each  set  con- 
sists of  two  differentiable  subsets  labelled  X  and  Y.    The  format  of  the 
outputs  is  as  described  earlier,  with  points  being  output  as  X,  0,  8,  E, 
or  J.     The  vector  C  in  the  discriminant  functions  is  the  centering  vector 
defined  in  Section  III.    Coordinate  axes  are  drawn  in.    The  points  picked 

for  iteration  are  circled,  and  the  direction  for  the  iteration  is  indicated 
i 

by  an  arrow  extending  from  the  circle. 

The  tuples  of  numbers  which  make  up  each  subset  were  generated 
on  the  IBM  360  as  uniformly  distributed  random  numbers  over  appropriate 
ranges.    In  each  case  the  number  of  elements  per  subset  depended  on 
the  dimensionality  of  the  problem;  for  n  dimension  each  subset  should 
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have  at  least  2     members.    This  lower  bound  is  set  to  insure  some 
reasonable  degree  of  definition  for  the  subsets  in  the  n  space.    The 
particulars  of  each  data  set  and  the  results  of  the  sample  runs  will  now 
be  discussed  sequentially. 

CASE  1 .     The  first  problem  to  which  the  algorithm  was  applied  in- 
volved the  separation  of  two  three  dimensional  unit  cubes.    Both  the  X 
and  Y  subsets  consisted  of  thirty-three  tuples.     If  X=(x    ,x_  ,x„)  then  x. 
for  1=1,3  was  generated  as  a  random  number,  uniform  (0,1);  if 
Y=(y,  ,V9  ,yJ  the  y.  for  i=l  ,3  was  generated  uniform  (1,2).    Thus  the  X 
'subset  can  be  thought  of  as  being  drawn  from  the  unit  cube  with  one 
vertex  at  the  origin  and  lying  in  the  three-space.     The  Y  subset  was  drawn 
from  the  unit  cube  displaced  diagonally  from  the  origin  by  the  X  cube. 
The  two  cubes  have  the  vertex  (1,1,1)  as  a  point  in  common;  however  in 
accordance  with  the  requirement  that  the  sets  be  disjoint,  (1,1,1)  was 
considered  to  be  a  possible  Y  value. 

Figures  VI- 1-1  through  VI- 1-3  show  a  three  step  sequence  in  the 
application  of  the  algorithm.    In  Figures  VI- 1-1  and  VI- 1-2  the  Y 
(indicated  as  O's)  and  X  subsets  are  mixed;  however  satisfactory  separation 
is  obtained  in  Figure  VI-1-3.    The  values  of  the  A  and  B  vectors  are  as 
follows: 

Figure  VI- 1-1        A(l)=-.5507  B(l)=-1.1837 

A(2)=    .5069  B(2)=  -.7223 

A(3)=-.1261  B(3)=   1.4448 
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Figure  VI- 1-1.     CASE  1 ,  FIRST  STEP 
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Figure  VI- 1-2.     CASE  1,  SECOND  STEP 
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Figure  VI- 1-3.     CASE  1,  THIRD  STEP 
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Figure  VI- 1-4.     PHOTO  OF  CASE  1,  FIRST  STEP 
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Figure  VI- 1-5.    PHOTO  OF  CASE  1,  SECOND  STEP 
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Figure  VI- 1-6.     PHOTO  OF  CASE  1,  THIRD  STEP 
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Figure  VI- 1-2        A(l)=-.3283  B(l)=-1.1837 

A(2)=    .5668  B(2)=   -.7223 

A(3)=-.0842  B(3)=    1.4448 

Figure  VI- 1-3        A(l)=    .5713  B(l)=-1.1837 

A(2)=    .5229  B(2)=   -.7229 

A(3)=    .6769  B(3)=    1.4448 

From  Figure  VI- 1-3  a  candidate  discriminant  function  can  be  clearly- 
seen  as 

f(Z)  =  A-  (Z-C)/TMAX 
where  A  has  the  component  values  listed  for  Figure  VI- 1-3  above,  and 
TMAX=1.519  8.     Thus  if  f(Z)  >    0,  Z  is  classified  as  a  Y;  if  f(Z)  <  0,  Z  is 
classified  as  an  X. 

Figures  VI- 1-4  through  VI-1-6  are  provided  to  show  the  degree  of 
correlation  between  the  view  on  the  AGT  screen  and  the  resulting  printed 
output.     Figures  VI- 1-4,  VI- 1-5,  and  VI-1-6  correspond  respectively  to 
the  outputs  designated  Figures  VI- 1-1,  VI-1-2,  and  VI-1-3;  the  high 
degree  of  similarity  between  the  related  figures  is  clear.    The  results 
of  this  comparison  can  be  readily  extrapolated  to  the  remaining  problems 
of  this  section  and  to  applications  of  the  algorithm  in  general;  therefore, 
photographic  evidence  for  the  remaining  cases  need  not  be  included. 

CASE  2  .    For  the  second  application  the  dimensionality  of  the  problem 
increased  from  three  to  six.    The  Y  and  X  subsets  were  composed  of 
seventy-six  tuples  generated  from  unit  cubes  oriented  in  six-space  in 
the  same  fashion  as  the  unit  cubes  of  case  one  were  in  three-space. 
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Thus  a  sample  from  the  Y  subset  would  be  Y=(y.  ,y„  ,y0  ,y.  ,y_  ,yJ  where 

1      Z       o      4       b       b 

l*y£2  for  i=l,6;  and  for  X  the  correspondence  would  be  X=(x    ,x    ,x    ,x    , 

xr,xj,   O^x.^1  for  i=l,6.     Figures  VI-2-1  through  VI-2-3  depict  the 
5      b  l 

separation  algorithm  applied  to  this  problem.     Note  that  except  for  the 
greater  number  of  data  points  involved,  the  complexity  of  the  problem  is 
not  increased  from  the  three  dimensional  problem  of  case  one.    The  A 
and  B  vectors  for  these  three  views  are  as  follows 


Figure  VI-2-1 

A(l)=   1.3448 

B(l)= 

.0868 

A(2)=      .3701 

B(2)=- 

-1.4109 

A(3)=-1.9594 

B(3)= 

-.7244 

A(4)=      .7525 

B(4)=- 

-1.0185 

A(5)=    1.3047 

B(5)= 

2.9777 

A(6)=   -.7075 

B(6)= 

.4225 

Figure  VI-2-2 

A(l)=    1.3448 

B(l)= 

.6450 

A(2)=      .3701 

B(2)= 

-.5899 

A(3)=-1.9594 

B(3)= 

.0692 

A(4)=      .7525 

B(4)= 

-.3409 

A(5)=    1.3047 

B(5)= 

3.0208 

A(6)=  -.7075 

B(6)= 

.7001 

Figure  VT-2-3 

A(l)=    1.3448 

B(l)= 

.7298 

A(2)=      .3701 

B(2)= 

.4801 

A(3)=-1.9594 

B(3)= 

.5716 

A(4)=      .7525 

B(4)= 

.5345 

A(5)=    1.3047 

B(5)= 

3.0897 

A(6)=  -.7075 

B(6)= 

.7335 
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Figure  VI-2-1.     CASE  2,  FIRST  STEP 
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Figure  VI-2-2.     CASE  2,  SECOND  STEP 
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Figure  VI-2-3.     CASE  2,  THIRD  STEP 
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From  Figure  VI-2-3  a  reasonable  discriminant  function  would  be 

f(Z)=B'(Z-C)/SMAX 
where  B  has  component  values  listed  for  Figure  VI-2-3  above  and 
SMAX=3  .  1940.     Thus  a  vector  Z  for  which  f(Z)  >  0  would  be  assigned  as 
a  Y;  if  f(Z)  <  0  the  result  would  be  classification  of  the  Z  as  an  X. 

CASE  3  .    In  the  last  problem,  separation  was  obtained  between 

two  sets  of  six-dimensional  data.    The  aim  of  the  problem  in  case  three 

is  to  complicate  the  situation  by  generating  two  six-dimensional  subsets, 

labelled  Y  and  X,  each  of  which  can  further  be  described  as  consisting 

of  two  subsets.    The  purpose  of  the  application  is  to  see  if  the  algorithm 

would  cause  a  separation  to  appear  not  only  between  the  Y  and  X  sets, 

but  also  within  these  sets.    In  a  sense  this  problem  is  a  combination  of 

a  pattern  classification  application  and  a  cluster  analysis  problem.    The 

Y  and  X  sets  each  consist  of  one  hundred  and  forty-six  tuples;  these  sets 

were  additionally  generated  as  two  subsets  of  seventy  members  each. 

Thus  the  first  seventy  points  in  the  Y  subset  were  generated  as  Y=(y    ,y~  , 

y0  ,y.  ,y_  ,y_)  where  0  ^  y.  <  I  for  1=1,4  and  1  -  y.  —  2  for  i=5,6 .    In 
£      4       b       b  l  l 

the  second  subset  of  Y's,  0  —  y.  —  1  for.  i- 1.6.    The  first  seventy  X's  were 

typically  X=(x,  ,x.  ,x_  ,x  .  ,xr  ,x.)  where  0  ^  x.  <  1  for  i=l ,  5  and 
1      I      o      4      5      b  l 

1  ^  x_  ^  2;  for  the  second  X  subset,  0  ^x.   <  1  for  i=l,4,   1   ^  xr  ^  2, 
b  lb 

and  0  ^  x_  ^    1.    As  in  case  two  the  increased  complexity  of  the  situation 
b 

is  reflected  immediately  by  the  number  of  points  involved;  the  computa- 
tional aspects  are  unaffected.    Once  again  only  three  figures  are  required 
to  display  separation.    The  A  and  B  vectors  assume  values  as  indicated 
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below: 


Figure  VI-3-1 

A(l)=      .0503 

B(l)=      .1174 

A(2)=   -.3371 

B(2)=   -.5137 

A(3)=      .5022 

B(3)=    1.4294 

A(4)=    1.1334 

B(4)=   -.5413 

A(5)=      .6375 

B(5)=      .6996 

A(6)=    1.5314 

B(6)=-1.2121 

Figure  VI-3-2 

A(l)=    1.6166 

B(l)=      .1174 

A(2)=-1.7057 

B(2)=  -.5137 

A(3)=      .6496 

B(3)=    1.4294 

A(4)=   1.2731 

B(4)=  -.5413 

A(5)=      .9450 

B(5)=      .6996 

A(6)=    1.9902 

B(6)=-1.2121 

Figure  VI-3-3 

A(l)=    1.6166 

B(l)=-1.5022 

A(2)=-1.7057 

B(2)=-2.1664 

A(3)=      .6496 

B(3)=    1.4117 

• 

A(4)=    1.2731 

B(4)=   -.3975 

A(5)=      .9450 

B(5)=      .6192 

A(6)=    1.9902 

B(6)=-1.4552 

The  discriminant  function  will  obviously  be  of  a  more  complicated  nature, 
judging  from  the  appearance  of  Figure  VI-3-3.     The  following  function 
appears  to  be  satisfactory: 

f(Z)=(A'  (Z-O/TMAX+B  -  (Z-C)/SMAX)(A  •  (Z-C)/TMAX-B  •  (Z-O/SMAX) 
where  A  and  B  assume  values  indicated  for  Figure  VI-3-3,  TMAX=7.3396, 
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Figure  VI-3-1.     CASE  3,  FIRST  STEP 
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Figure  VI-3-2.     CASE  3,  SECOND  STEP 
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Figure  VI-3-3.     CASE  3,  THIRD  STEP 
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and  SMAX=7.5366.    For  this  problem,  f(Z)  >0  would  result  in  Z  being 
classified  as  an  X;  f(Z)  <  0  naturally  would  result  in  Y  classification. 

CASE  4 .    In  the  first  data  sets  of  this  section  the  points  were 
generated  from  multi-dimensional  cubes.    In  case  four  a  different  orienta- 
tion is  used,  this  time  in  four  dimensions.     Twenty-four  tuples  are 

generated  for  each  of  the  X  and  Y  sets .    The  Y  set  generated  so  that  a 

2  2  2  2    >> 

typical  point  Y=(y    ,y  ^  'Y3  'Y4)  satisfies  2  £  y    +  y  ^  +  y^  +  y^  ^  4 . 

2         2         2         2 
Similarly  an  X  point  satisfies  x     +  x_  +  x„  +  x     -  1.     Thus  the  geometric 

"picture"  is  that  of  a  sphere  surrounded  by  a  shell.    For  this  application, 

two  sets  of  figures  are  provided  to  emphasize  the  importance  of  the  initial 

values  assigned  to  the  A  and  B  vectors.    Figures  VI-4-1  through  VI-4-5 

display  the  results  of  making  an  unfortunate  choice  of  initial  values  for 

A  and  B.    Thirty  iterations  were  carried  out  from  this  initial  choice; 

however,  no  significant  improvement  on  the  situation  as  indicated  by 

Figure  VI-4-5  was  made.    The  values  of  the  A  and  B  vectors  were  as  follows: 


Figure  VI-4-1 

A(l)= 

.5780 

B(l)=-    .5000 

A(2)= 

.5780 

B(2)=-    .5000 

A(3)= 

.5780 

B(3)=    1.0000 

A(4)= 

.0000 

B(4)=    1.0000 

Figure  VT-4-2 

A(l)= 

.6680 

B(l)=-    .5000 

• 

A(2)=- 

.2623 

B(2)=-    .5000 

A(3)=- 

.0208 

B(3.)=    1.0000 

A(4)= 

.0900 

B(4)=    1.0000 
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Figure  VI-4-3 

A(l 

)=   .6680 

b(i; 

l=-  .0280 

A(2, 

!=-  .2623 

B(2, 

1=   .0917 

• 

A(3 

>=-  .0208 

B(3] 

1=   .4002 

A(4, 

1=   .0900 

B(4] 

1=   .7390 

Figure  VI-4-4 

A(l 

)=   .6480 

b(i: 

l=-  .0280 

A(2, 

)=   .1963 

B(2] 

1=   .0917 

A(3, 

>=-  .7196 

B(3, 

)=   .4002 

A(4, 

1=   .6110 

B(4, 

1=   .7390 

Figure  VI-4-5 

A(i: 

>=   .6680 

B(l 

►=-  .0280 

a(2; 

>=-  .2623 

B(2] 

1=   .0917 

a(3; 

l=-  .0208 

B(3, 

1=   .4002 

a(4; 

1=   .0900 

B(4] 

)=   .7390 

Rather  than  attempt  to  find  a  satisfactory  discriminant  function  for  the 
situation  depicted  in  Figure  VI-4-5,  the  approach  is  to  specify  new  A  and 
B  vectors  and  to  start  over.     In  fact,   starting  with  A=(l ,  1 , 1 , 1)  and 
B=(l  ,-1 , 1  ,-1) ,   separation  appeared  within  twenty  iterations.     Figures 
VI-4-6  through  VE-4-10  depict  the  situation  after  the  first,  fifth,  tenth, 
fifteenth,  and  twentieth  iterations  respectively.    The  A  and  B  vectors  for 
the  figures  are: 

Figure  VI-4-6       A(l)=   1.0000  B(l)=   1.0000 

A(2)=    1.0000  B(2)=-1.0000 

A(3)=    1.0000  B(3)=    1.0000 

A(4)=    1.0000  B(4)=-1.0000 
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Figure  VI-4-1.     CASE  4,  FIRST  STEP 
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Figure  VI-4-2  .     CASE  4,  SECOND  STEP 


52 


**«***«»#*«***********#«*-»«***< 

}#*********#»**##*»**-l(-**##*##*# 

* 

* 

* 

♦ 

# 

* 

* 

♦ 

*                                   0 

0 

♦ 

*                                                                     0 

♦ 

* 

* 

0 

0 

♦ 

* 

* 

* 

♦ 

« 

; 

♦ 

» 

XX 

« 

♦  0 

♦ 

* 

X 

X 

X 

♦ 

* 

X 

E 

♦ 

*  0                            0 

♦  000                                                   (yV 

* 

w 

J                                                                                                                                                                  w 

X 

U 

*                                 0 

XE 

* 

* 

X 

X 

« 

■ 

X 

X 

« 

* 

X 

* 

« 

X        x 

* 

*                                                       0 

« 

# 

« 

*                                               0 

♦ 

# 

« 

# 

* 

« 

# 

« 

# 

* 

# 

*                                             0 

- 

♦ 

« 

« 

»                                             .8         0 

# 

# 

0 

* 

*#•*»#♦*****#♦  +  *»♦#*»#  » .**»***««§  #♦»###*#*#*»♦»#<(.###*♦♦*.*#♦«♦*♦ 


Figure  VI-4-3.     CASE  4,  THIRD  STEP 
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Figure  VI-4-4.     CASE  4,  FOURTH  STEP 
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Figure  VI-4-5.     CASE  4,  FIFTH  STEP 
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Figure  VI-4- 

-7 

A(l)=   2.1930 

B(l)=      .9800 

A(2)=   2.1246 

B(2)=-    .5413 

A(3)=-    .8528 

B(3)=      .3012 

A(4)=-    .2470 

B(4)=-    .4790 

Figure  VI-4- 

-8 

A(l)=    1.8910 

B(l)=    1.5190 

A(2)=      .4019 

B(2)=-    .4473 

A(3)=-    .5144 

B(3)=    1.7293 

A(4)=      .5510 

B(4)=-    .4620 

Figure  VI-4- 

-9 

A(l)=    1.5890 

B(l)=      .8490 

A(2)=-1.3207 

B(2)=-    .9287 

A(3)=-    .1760 

B(3)=   3.3314 

A(4)=    1.3490 

B(4)=-    .1040 

Figure  VI-4- 

-10 

A(l)=   2.1580 

B(l)=    1.6860 

A(2)=-    .3294 

B(2)=-1.0144 

A(3)=-    .6366 

B(3)=   2.3048 

A(4)=    1.4610 

B(4)=-1.2420 

The  situation  in  Figure  VI-4-10  provides  another  opportunity  for  displaying 
the  flexibility  possible  in  choosing  a  discriminant  function.     The  function 
in  this  case  is  circular  with  the  equation: 

f(Z)=  (A-  (Z-O/TMAX-1/3)2  +  (B-(Z-C)/SMAX)2-l/4 
where  A  and  B  have  the  values  specified  for  Figure  VI-4-10,  TMAX=3.3577, 
and  SMAX=5  .  8823  .     Clearly  if  f(Z)  <  0  ,  Z  is  classified  as  an  X;  if  f(Z)  >  0  , 
Z  is  classified  as  a  Y. 
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Figure  VI-4-6.     CASE  4,  SIXTH  STEP 
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Figure  VI-4-7.     CASE  4,  SEVENTH  STEP 
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Figure  VI-4-8.     CASE  4,  EIGHTH  STEP 
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Figure  VI-4-9.     CASE  4,  NINTH  STEP 
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Figure  VI-4-10.     CASE  4,  TENTH  STEP 
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It  should  be  noted  that  the  separation  found  was  a  result  of  the 
particular  samples  used.    In  reality  the  Y  elements  must  have  formed  a 
ring  instead  of  a  shell.     Clearly  there  is  no  possible  orientation  for  the 
A-B  plane  into  which  the  sphere  and  shell  could  be  projected  so  as  to 
display  separation.    This  is  a  case  where  the  samples  were  in  fact  too 
small  to  allow  a  discriminant  function  to  be  derived  which  would  be 
applicable  to  the  parent  populations. 

Cases  five  and  six  represent  applications  of  the  separation  program 
to  "real  world"  data  sets.    The  complexities  found  in  actual  data  are  re- 
flected by  the  apparent  reduction  in  the  performance  level  of  the  algorithm 

CASE  5  .    The  data  set  for  this  cluster  analysis  problem  is  the 
Jamaican  "Fossil"  data  mentioned  in  the  discussion  of  Chemoff's  method 
in  section  two  of  this  paper.    The  data  set,  in  the  form  of  eighty-two 
eight-tuples,  was  drawn  from  Wright  &  Switzer's  paper  "Numerical 
Classification  Applied  to  Certain  Jamaican  Eocene  Nummulitids"    6     . 

As  indicated  earlier  in  this  paper,  the  data  was  represented  as  a 
single  set  of  points  Y   (points  on  the  output  are  indicated  only  as  O's  or 
8's),  and  iterations  were  made  based  on  apparent  separation;  the  aim  is 
to  emphasize  this  separation.     Twelve  iterations  were  made  starting  with 
A=(l,  1,1, 1,1, 1,1,1)  and  B=(l,-1, 1,-1, 1,-1, 1,-1).    Figures  VI-5-1 
through  VI-5-6  were  output  on  the  first,   second,  third,  fourth,  ninth, 
and  twelfth  iterations  respectively.    Since  no  discriminant  function  is 
generated  in  a  cluster  analysis,  the  values  of  the  A  and  B  vectors  have 
no  final  importance.    On  the  other  hand,  the  precise  identity  of  each 
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Figure  VI-5-1.     CASE  5,  FIRST  STEP 
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Figure  VI-5-2.     CASE  5,  SECOND  STEP 
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Figure  VI-5-3.     CASE  5,  THIRD  STEP 
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Figure  VI- 5-4.     CASE  5,  FOURTH  STEP 
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Figure  VI- 5-5.     CASE  5,  FIFTH  STEP 
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Figure  VI- 5-6.     CASE  5,  SIXTH  STEP 
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Y  point  in  the  output  is  now  important;  this  information  is  deducible  from 
the  point  coordinates  which  are  also  part  of  the  program  output.     The 
author  found  that  the  graphic  output  provided  indications  generally  of  two 
clusters;  one  consisted  of  points  one  through  thirty-eight  and  seventy- 
three  through  eighty-two,  and  the  other,  of  points  thirty-nine  through 
seventy-two.    These  clusters  are  indicated  in  Figure  VI- 5-5,  which 
represents  the  situation  after  nine  iterations.     Comparing  these  results 
to  those  of  Wright  &  Switzer  who  found  three  clusters,  it  appears  that 
the  present  method  failed  to  differentiate  between  their  groups  one  and 
three.     The  clusters  found  by  Wright  &  Switzer  have  been  indicated  in 
Figure  VI-5-6  (which  represents  the  situation  after  three  more  iterations), 
identifying  them  by  the  coordinates  of  their  respective  members  rather 
than  by  operator  ovservation.    Inspection  of  the  positions  of  these  groupings 
indicates  that  the  clustering  scheme  proposed  by  Wright  &  Switzer  might 
eventually  be  attained  or  at  least  approached  by  the  algorithm's  appli- 
cation.   An  important  point  is  that  Wright  &  Switzer  required  ninety-six 
iterations  to  arrive  at  their  solution;  all  results  of  the  present  method 
were  found  in  less  than  thirteen  iterations. 

CASE  6 .    The  second  "real  world"  problem  to  be  approached  with 
this  algorithm  is  a  thirteen  dimensional  pattern  classification  problem. 
Each  thirteen-tuple  represents  a  set  of  medical  readings  taken  from 
coronary  patients  at  Letterman  Hospital  in  the  case  of  the  Y's,  and  taken 
from  personnel  receiving  physical  examinations  at  Fort  Ord  Hospital  for 
the  X's.    The  Letterman  patients  were  known  to  have  been  suffering 
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from  coronary  disease.     No  such  concrete  statement  can  be  made  for  the 
Fort  Ord  examinees.    For  the  latter  it  was  not  absolutely  known  that  they 
did  not  have  heart  disease.     This  separation  along  with  the  appearance  of 
the  data  indicated  that  if  separation  existed,  it  would  probably  not  be 
well  defined.     Certain  of  the  features  in  each  thirteen-tuple  represented 
measurable  quantities  such  as  age  or  weight.    Other  features  such  as 
race,  medical  history,  and  smoking  habits  were  classifiers  rather  than 
measurements . 

The  application  of  the  algorithm  failed  to  achieve  a  clear  cut 
separation  between  the  two  sets;  however  as  the  sequence  of  Figures 
VI-6-1  through  VI-6-6  indicates,   relatively  high  concentrations  of  Y's 
and  X's  were  attained.    The  A  and  B  vectors  for  Figure  VI-6-6  are  as 
follows: 


Figure  VT-6-6 


A(l)=   .0000 

B(l)=-  .9273 

A(2)=37.0000 

B(2)=-12.4091 

A(3)=-20.0000 

B(3)=  7.0818 

A(4)=   .0000 

B(4)=38.3364 

A(5)=-1.0000 

B(5)=   .6455 

A(6)=-2.0000 

B(6)=-  .69  09 

A(7)=-1.0000 

B(7)=   .1636 

A(8)=-2.0000 

B(8)=-1.4636 

A(9)=   .0000 

B(9)=   .6636 
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Figure  VI-6-1.     CASE  6,  FIRST  STEP 
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Figure  VI-6-2  .     CASE  6,  SECOND  STEP 
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Figure  VI-6-3.     CASE  6,  THIRD  STEP 
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Figure  VI-6-4.     CASE  6,  FOURTH  STEP 
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Figure  VI-6-5.     CASE  6,  FIFTH  STEP 
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Figure  VI-6-6.     CASE  6,  SIXTH  STEP 
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A(10)=        .0000  B(10)=        .2545 

A(lD=-    1.0000  B(ll)=        .5818 

A(12)=-10.0000  B(12)=   69.0818 

A(13)=      3.0000  B(13)=-    8.4455 

Obviously  no  clear  cut  recommendation  for  a  discriminant  function  is 
justified  by  Figure  VI-6-6;  however,  the  X  data  can  be  seen  to  be  concen- 
trated in  the  third  quadrant  of  the  figure.    This  information  should  be 
useful  in  making  tentative  discriminant  judgements.    A  possible  discrim- 
inant function  would  be: 

f(Z)=B  ■  (Z-C)/SMAX  +  (A-  (Z-C)/TMAX)/10    for  A*  (Z-C)  <  0 
=B-(Z-C)/SMAX  +  3(A-(Z-C)/TMAX)/4    for  A-(Z-C)  >  0 
This  function  classifies  thirty-seven  out  of  fifty  X  points  and  fifty-three 
out  of  sixty  Y's  correctly.     This  result  is  not  satisfactory  as  the  results 
obtained  in  the  first  four  cases;  however  it  is  probably  more  reasonable 
to  use  it  than  to  conduct  the  great  number  of  iterations  likely  to  be  required 
to  find  separation  (if  it  exists)  and  a  more  suitable  discriminant  function. 
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VII.     CONCLUSIONS  AND  RECOMMENDATIONS 

Based  on  the  results  noted  using  the  generated  data  sets,  it  can  be 
concluded  that  the  separation  algorithm  proposed  by  Professor  Cover  is 
satisfactorily  implemented  by  the  program  written  for  this  paper.     Provided 
that  there  exists  an  orientation  for  the  A-B  plane  in  which  the  projections 
of  the  data  points  are  separated,  the  algorithm  furnishes  an  avenue  for 
describing  this  orientation.     In  both  real  world  applications  such  an 
orientation  was  approached,  but  not  actually  found.    It  is  apparent  that 
dimensionality  per  se  is  not  a  problem  for  the  separation  algorithm.     The 
display  for  the  thirteen-dimensional  hospital  problem  was  no  more  complex 
than  that  for  the  generated  set  of  four  six-dimensional  cubes.    Dimension- 
ality can  complicate  implementation  of  the  algorithm  when  the  range  of 
satisfactory  A-B  plane  orientations  is  small  relative  to  the  sample  space. 
A  small  range  requires  in  general  a  greater  number  of  iterations.     Once 
separation  is  achieved,  definition  of  the  discriminant  function  is  quite 
simple.    Using  the  printed  output  of  the  program,  the  discriminant  function 
can  be  drawn  in  and  its  equation  found  as  a  polynomial  curve. 

An  important  feature  of  this  separation  algorithm  is  that  it  makes 

no  distributional  assumptions  or  requirements  on  the  data.     This  is 

valuable  indeed  when  working  with  "real  world"  situations.    Additionally 

the  chain  of  operations  on  a  data  point  are  simple  and  easy  to  follow. 

Once  a  discriminant  function  has  been  defined,  it  is  easily  applied  to 

new  data  points . 
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As  a  final  comment,  the  importance  of  the  man-machine  approach 
of  this  algorithm  cannot  be  overestimated.    Increased  familiarity  with  the 
data  seems  to  reduce  the  number  of  iterations  required  for  the  achievement 
of  separation.    As  the  same  data  set  was  used  repeatedly,  it  was  obvious 
that  the  operator  got  a  "feel"  for  the  movement  observed.    Also  the  man's 
importance  is  emphasized  when  it  is  realized  that  he  makes  all  decisions 
and  judgements  during  the  application.    It  was  found  that  the  choice  of 
X  and  Y  as  the  identifiers  for  the  two  classes  in  the  problem  was  unfor- 
tunate.   The  similarity  in  appearance  of  these  two  letters  caused  some 
confusion  in  the  operator.    It  was  to  alleviate  some  of  this  confusion 
that  the  labelling  scheme  for  the  output  used  X's  and  O's  instead  of 
X's  and  Y's. 

Recommendations  for  future  work  with  the  separation  algorithm  are 
aimed  at  allowing  the  operator  to  gain  a  better  grasp  of  the  situation  he 
views  on  the  screen  of  the  AGT.     First,  the  video  display  should  be 
changed  from  the  present  X-Y  scheme  to  some  other  less  confusing  scheme 
such  as  an  X-O  display. 

Second,  a  feature  should  be  available  in  the  program  to  indicate 
when  specified  densities  of  points  are  shown  on  the  screen.     For  example, 
if  when  the  density  of  X's  in  a  part  of    the  screen  reaches  a  certain  level, 
the  video  of  these  X's  either  goes  out  of  focus  or  assumes  a  differentiable 
intensity. 
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Third,  the  application  of  the  algorithm  might  be  extended  by  making 
the  projection  onto  a  surface  other  than  a  hyperplane.     This  modification 
requires  significant  alteration  to  the  work  presented  here. 

Fourth,  it  might  be  helpful  to  the  operator  to  have  video  in  the 
shape  of  an  arrow  appear  by  the  point  to  be  moved,  indicating  direction 
of  movement.    Also  since  each  member  of  the  data  set  has  an  identifying 
number,  it  might  be  helpful  when  choosing  a  point  to  have  the  number  of 
that  point  appear  on  the  screen. 
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APPENDIX  I   PRGGRAN  LISTING 


SET  UP  CF  ARRAYS,  INITIALIZING 
CONSTANTS,  AND  DATA  INPUT 


INTEG 

CIMEN 

CIKEN 

CINEN 

DIMEN 

DIMEN 

DIMEN 

DIMEN 

DATA 

NAf'EL 

IDEV= 

CLTPU 

INPUT 

CALL 

I  F  (  IE 

TT  =  0. 

N  =  50 

N  =  6C 

K=13 

READ( 

IF(M. 

READ( 

FCPMA 


ER  ALFH(8) ,PIC(61,37) 

SIGN  ITDIR(3) 

SIGN  IGDIR(3) 

SIGN  T0T(13 )  ,0AV<  13) 

SIGN  IMAGO01)  ,IFRAM(30L) 

SIGN  JMAG(301) , JFRAM(301) 

SICN  X(50,13),Y(6C,13),XX1(5C) ,XX2 (50) 

SIGN  YYK60)  ,  YY2(  60) 

8)/  1H0, 1H8, 1HX,  IFE,1HJ  , 

IER 


A(  13)  ,B(  13) 


YYK60) 
(ALPH(  I)  ,1  =  1 
1ST  A,B,IDEV 
1 

T(101)  'TYPE 
(101) 
DTINITt IDEV, 


IF  ,  1H*,1HC/ 


R.NE 
01 


IDEV=2*  AND  C/R  IF  USING  AGT  2' 
ITDIR,3, IER) 


0)  OUTPUT(lOl)  IER,'DTINIT« 


5,2187)((Y(I,J),J=1,K),I=1,N) 
EQ.O)  GO  TO  167 
5f2l87)((X(IfJ),J=lfK)iI=lff) 
T(F4.1,12F5 .0) 


CENTRALIZING  DATA 

DC  170  J=1,K 
TCT(J)=0.0 

CCNTINUE 

DC  171  11=1, N 

DC  172  JJ=i,K 

TCT( JJ)=TCT(JJ)+Y( I  I, JJ) 

CONTINUE 

CONTINUE 

DC  173  II=l,M 

DC  174  JJ=1,K 

TCT(JJ)=TCT(JJ)+X( II, JJ) 

CCNTINUE 

CONTINUE 

DEN=FL0AT(M+N) 

DC  175  JA=1,K 

DAV(JA)=T0T( JA)/DEN 

WPITE(6  ,181)  JA, DAV(JA) 

F0RMAT( 5X, 'DATA  CONP  .  AVE 

CCNTINUE 

N 


NC 


13 


,F10.3) 


IK)-DAV( IK) 


,K 


DC  176  IJ=1 
DO  177  IK=1,K 
Y(  IJ,IK)=Y(  IJ, 

CCNTINUE 

CCNTINUE 

CC  178  IJ=1 

DC  179  IK=1 

X( IJ,IK)=X( IJ,IK)-DAV( IK) 

CCNTINUE 

CCNTINUE 

WRITE(6,  190) 

FCPMAT(5X, 'CENTRALIZED  DATA',/) 

DO  191  IJ=1,N 

DC  192  IK=1,K 

WPITE(6,193) 

FORMAT*  5X, »Y( ' , 13 

CCNTINUE 

hPITE{6,194) 

FCRMAT( 5X,/) 

CCNTINUE 


IJ,  IK,Y(IJ,IK) 


',  ', 


13,' )=« ,F10.3) 
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c 
c 
c 


197 
196 

195 

200 


450 
30 


10 


C 

c 
c 


15 
12 


17 
14 

182 


C 

c 


11 


DC  195  IJ=1,M 

DO  196  IK=1,K 

WRITE(6,197)  IJ,  IK,X(  IJ, IK) 

FCRMAT(5X,'X( ',I3,',',I3,')=',F10.3) 

CONTINUE 

WPITE(6 ,194) 

CCNTINUE 

WRITE(6,20C) 

FCPMAT(5X,///) 

INPUT  COMPONENT  VALUES  FOR  THE  A  AND  B  VECTORS 

CALL  GINPUT( IDEV, ITCIR,1) 

CC  30  JJ=1,K 

WPITE(6,450)  JJ, A( JJ)  ,JJ,B( JJ) 

F0RMAT(5X, 'A(',I3,')=',F6.3,'  ;B(',I3,')=',F6.3) 

CCNTINUE 

KST  =  0 

KRT=0 

CALL  DGINIT  (  I  DEV  ,  IGC  IR  ,  3  ,  I  ER  ) 

IF(IER.NE.O)  OUTPUT  (101)  IER, 'DGINIT' 

CALCULATION  OF  DATA  POINT  PROJECTIONS 


SMAX=0 

TPAX=0 

SLMYl= 

SUMY2= 

DC  12 

YYKI) 

YY2(I) 

DC  15 

YYK  I) 

YY2(  I) 

CONTIN 

SLMYl= 

SLPY2= 

CCNTIN 

SUMX1= 

SLMX2= 

DC  14 

XXKL) 

XX2(L) 

DC  17 

XXKL) 

XX2(L) 

CONTIN 

SUMX1= 

SUPX2= 

CCNTIN 

AV1=(S 

AV2=(S 

WRITE( 

FORMAT 


*Y(  I 
*Y(  I 


»J) 

tJ) 


I  ) 
I  ) 


.0 

.0 

0.0 

0.0 

1=1, N 

=  C0 

=  0.0 

J  =  1,K 

=YY1(I )+A( 

=YY2(I )+B( 

UE 

SUMYl  +  YYK 

SUMY2+YY2( 

UE 

0.0 

0.0 

L  =  1,M 

=  0.0 

=  0.0 

J  =  l,K 

=XXl(L)+A( 

=XX2(L)+B( 

UE 

SUMX1  +  XXHL) 

SLMX2+XX2( L) 

UE 

UMY1+SUMX1 ) /FLOAT (M+N) 

Ut*Y2  +  SUNX2)  /FLCAT(M+N) 

6,182)  AV1.AV2 

(5X,«SCAL.  AVE  1=',F10.6,'   SCAL.  AVE  2=',F10.6) 


J)*X(L,J) 
J)*X(Lf J) 


SCALING  THE  PROJECTIONS 

DC  11  1=1, N 
YYK  I)  =  YYl(I  J-AVi 
YY2(  I)=YY2( I )-AV2 
GA  =  ABS(  YYK  I  )  ) 
GE=ABS(YY2( I ) ) 
IF(GA.GT.TMAX)  TM£X=GA 
IF(GB.GT.SMAX)  SMAX=GB 
CCNTINUE 
DC  13  L=1,M 
XXKL)=XX1(L)-AV1 
XX2(L)=XX2(L )-AV2 
2A=ABS(XX1 (L) ) 
ZE=ABS( XX2(L) ) 
IF(ZA.GT.TMAX)  TMAX=ZA 
IF(ZB.GT.SMAX)  SMAX=ZB 
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47 


C 
C 

c 


13         CONTINUE 

KN=5*N 

KK=KN+1 

JL  =  0 

IFPAM( L  )=IHEAD(0, 101 

WPITE(6,6667)     TMAX,SMAX 
6667    F0RMAT(5X,  'TMAX='  ,F  L0.6,'        SMAX  =  »  ,F  L0.6  ) 

DRAWING    THE    Y    POINTS    ON    THE    AGT    SCREEN 

CC    47     1=2, KN, 5 

JL=JL+L 

FA=YYL(  JD/TMAX 

FB  =  YY2(  JD/SMAX 

WRITE(6,5001 )     JL,FA,JL,FB 

FC=FA+TT 

FD=FA-TT 

FE=FB+TT 

FF=FB-TT 

IFRAM(I  )=IPACK(FA, FB,0) 

IFRAM< I+1)=IPACK<FC,FE,1) 

IFRAM(I+2)=IPACK ( FA,FF,0) 

IFRAM(I+3)=IPACK(FA,FB,1) 

IFRAM( 1+4) =1  PACK (FD,FE,  L) 

CCNTINUE 

CALL  GRAPHO  ( I DE V  ,  I  FR AM , KK , 1 , IE P  ) 

IF(IER.NE.O)  OUTPUT  (L01)  IER.'GELKl1 

IN=6*M 

IF(  IN. ECO)  GO  TO  55 

JJ=6*M+1 

JN=0 

IMAG( 1) =IHEAD(0,10) 

DRAWING  THE  X  POINTS  ON  THE  AGT  SCREEN 

DO  45  1=2, IN, 6 

JN=JN+1 

TA=XX1( JN) /TMAX 

TB=XX2( JN)/SMAX 

WPITE(6  ,5000)  JN,TA,JN,TB 

TC=TA+TT 

TC=TA-TT 

TE=TB+TT 

TF=TB-TT 

IMAG( I )=IPACK(TA,TB,0) 

INAG( I+1)=IPACK(TC,TE,1) 

IMAG(I+2)=IPACK(TC,TF,0) 

IMAG(I+3)=IPACK(TD,TE, L) 

IMAG( I+4)=IPACK(TD,TF,0) 

IMAG( I+5)=IPACK(TA,T6, L) 

CCNTINUE 

CALL  GRAPHO  ( IDEV 


45 


C 
C 
C 


55 


C 

c 


IMAG, JJ,2,  IER) 
IF(IER.NE.O)  OUTPUT(IOL)  IER,'GeLK2' 

CHECKING  FOR  DEPRESSED  FUNCTION  SWITCHES 


CALL  FNS(  IDEVf ISW, IER) 
IF(IER.NE.O)  OUTPUT(IOL) 
IF(ISW.GE.O)  GO  TO  55 
IF(LLS(  ISW, 11 J.LT.O)  GO  TO  7 
IF(LLS(  ISW, 3) .LT.O)  GO  TO  57 
IF(LLS( ISW, 13).LT.O)  GO  TO  58 
IF(LLS( ISW, 14) .LT.O)  GO  TO  67 

,15). LT.O)  GO  TO  68 

1)  GO  TO  90 

1)  GO  TO  91 
16). LT.O)  GO 


IERj'GELKS' 


IF(LLS( 
IFCKRT. 
IF(KST. 
IF(LLS{ 
IF(LLS( 
IF(LLS( 


ISW 

EC 

EQ 

ISW 

ISW 

ISW 


GC  TO  55 


,17) 
,18) 


LT 
LT 


0) 
C) 


GO 
GO 


TO 
TO 
TO 


55 
55 
55 


CHOICE  CF  DIRECTION  CF  MOVE  WITH  Y  PICKED 


83 


90 


CALL  FNSC  IDEVflSW, IER) 


C 
C 

C 


91 


57 
L5L 


C 
C 
C 


601 
6CC 

58 
152 


C 

c 
c 


611 
610 

59 


IF( IER. NE.O) 
IF(  ISW.GE.O) 
IF(LLS(  ISW, 5) 
IF(LLS(  ISW, 6) 
IF(LLS( ISW, 7) 
IF(LLS(  ISW, 8)  .LT.O) 
GO  TO  5  5 


OUTPUT(lOl) 
GO  TO  55 
.LT.O)  GO 
.LT.O) 
LT.O) 


GO 
GO 
GO 


IERt'GELK?1 

TO  59 
TO  60 
TO  6  1 
TO  62 


CHOICE  OF  DIRECTION  OF  MOVE  WITH  X  PICKED 


CALL  FNS(  IDEV,  ISW,  IER) 
IF(IER.NE.O)  OUTPUT(lOl) 


IER, • GELK8' 


IF(  ISW.GE.O) 
IF(LLS( ISW, 5 
IF(LLS( ISW, 6 


GO  TO  55 
).LT.O)  GO 
) .LT.O) 


IFCLLSC ISW, 7) .LT.O) 


IF(LLS( ISW, 8 
GO  TO  5  5 

PICKING  A  Y 


)  .LT.O) 


GO 
GO 
GO 


TO  63 

TO  64 

TO  65 

TO  66 


IER  ) 
i 


151 


CALL  GRAPHR{ IDEV,  IFRAM,KK, 1, 

IF(IER.NE.O)  OUTPUT(lOl) 

IFCMODC  IGDIR(  1) , 8)  .EQ.O)  GO  TO 

CALL  GRAPHIC ICEV, JFRAM, 1, IER) 

IFCIER.NE.C)     0UTPUTC101)     IER,MGBLK' 

KEDS=0 

KST=0 

KRT  =  1 

KKK  =  0 

DC  600  1=2, KN, 5 

KKK=KKK+1 

IFCJFRAMC I ).NE. IFRAMC I ) )    GO 

IFC JFRAM C I +  1)  .NE. IFRAMC  1+1) 

.NE 

.NE 
NE 


TO  6C1 


IFCJFRAM  1+2) 
IFCJFRAMC 1+3) 
IFCJFRAM  1+4) 
GO  TO  600 
JEC=KKK 
GC  TO  55 
CONTINUE 

PICKING  AN  X 


IFRAMC 1+2) 
IFRAMC 1+3) 
IFRAMC 1+4} 


GO 
GO 
GO 
GO 


TO 
TC 
TO 
TG 


601 
601 
601 
601 


CALL  GR 
IFC  IER. 
I  F  ( MO  D  C 
CALL  GR 
IFC  IER. 
KEDS=0 
KST=i 
KRT  =  0 
LLL  =  0 
DC  610 
LLL=LLL 
IFCJMAG 
IFCJMAG 
IFCJMAG 
IFCJMAG 
IFCJMAG 
IFCJMAG 
GC  TO  6 
LED=LLL 
GC  TO  5 
CCNTINU 


IER) 
i 


APHRC ICEV,  IMAG, JJ, 2, 

NE.O)  0UTPUTC101) 

IGDIRC2) ,8) .EQ.O)  GO  TO  152 

APHIC IDEV, JMAG,2, IER) 

NE.O)  0UTPUTC101)  IER,'JGBLK' 


1  =  2 
+  1 
CI  ) 
CI  + 
CI  + 
C  1  + 
CI  + 
CI  + 
10 

c 
E 


IN, 6 


IMAGCD)  GO 

NE.  IMAGC  1  +  1) 

NE.  IMAGC 1+2) 

NE.  IMAGC 1  +  3) 

■  NE.  IMAGC 1+4) 


.NE 

1) 

2) 

3) 

4) 

5) .NE. IMAGCI+5) ) 


TO  611 

GO  TO 

GO 

GO 

GO 

GO 


TC 
TO 
TO 
TO 


611 
611 
611 
611 
611 


ITERATION  BY  MOVING  Y  TC  THE  RIGHT 

IFCJED.EC.O)  GO  TO  55 
WRITE(6,2000)  JED 
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2000  F0RMAT(5X,  'MOVING  RIGHT  USING  YCB,1!') 
DC  110  KJ=1,K 

A(KJ)=A(KJ)+Y( JED,KJ) 

WPITE(6,  1000)  KJ,A(KJ) 
110   CCNTINUE 

KEDS=1 

GC  TO  10 
C 

C      ITERATION  BY  MOVING  Y  TO  THE  LEFT 
C 

60  IF(JED.EC.O)  GO  TO  55 
WRITE(6,2001)  JED 

2001  F0PMAT(5X, 'MOVING  LEFT  USING  YfSIS,')') 
DC  111  KJ=1,K 

A(KJ)=A(KJ)-Y( JED,KJ) 

WRITE(6, 1000)  KJ»A(KJ ) 
ill   CCNTINUE 

KECS=2 

GC  TO  10 
C 

C      ITERATION  BY  MOVING  Y  UP 
C 

61  IF(JED.EQ.O)  GO  TO  55 
WRITE(6,2002)  JED 

2002  FCPMAT(5X,  'MOVING    UP    USING    YC  ,13,')') 
DO    112    KJ=1,K 

B(KJ)=B(KJ )+Y(JED,KJ) 
WRITE(6,1001 )     KJ,B(KJ) 

112  CONTINUE 
KECS=3 
GC    TO    10 

C 

C      ITERATICN  BY  MOVING  Y  DOWN 

C 

62  IF(JED.EQ.C)  GO  TO  55 
WRITE(6,2003)  JED 

2C03  FCRMATC 5Xf •MOVING  DOWN  USING  YCttI3»,Ji ) 
DO  113  KJ=i,K 
B(KJ)=B(KJ)-Y(JED,KJ) 
WRITE(6,l00l)  KJ,B(KJ) 

113  CCNTINUE 
KEDS=4 
GC  TO  10 

C 

C      ITEPATICN  BY  MOVING  X  TO  THE  RIGHT 

C 

63  IF(LED.EQ.O)  GO  TO  55 
WRITE(6,3000)  LED 

300C  FCRMAT(  5X,  'MOVING  RIGHT  USING  X(,,I3,,)M 
CC  120  KD=1,K 
A(KD)=A(KD)+X(LED,KC) 
WPITE(6» 10C0)  KD,A(KD) 

120  CONTINUE 
KECS=5 
GO  TO  10 

C 

C      ITERATICN  BY  MOVING  X  TO  THE  LEFT 

C 

64  IF(LEC.EG.O)  GO  TO  55 
Vn'RITE(6,3001  )  LED 

3001  FORMAT!  5X, 'MOVING  LEFT  USING  Xt't^,')1) 
DC  121  KD=1,K 
A(KD)=A(KD)-X(LED,KD) 
WRITE(6,1000)  KD,A(KD) 

121  CCNTINUE 
KEDS=6 
GC  TO  10 

C 

C      ITERATICN  BY  MOVING  X  UP 

C 

65  IF(LEC.EQ.O)  GO  TO  55 
WRITE(6,3C02)  LED 
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c 
c 
c 


2002 
122 

66 
3002 

122 
67 


827 

825 


704 


705 


7  1  A 


715 


777 


890 

888 

751 


FORMATC  5X,  'MOVING  LP  USING  XCn,')1) 

DC  122  KD=1,K 

B(KD)=B(KD)+X(LED,KD) 

WRITE(6,  1001)  KD,B(KD) 

CONTINUE 

KEDS=7 

GC  TO  10 

ITERATION  BY  MOVING  X  DCWN 

IF(LEC.EC.O)  GO  TO  55 

WRITE(6,3003)  LED 

FOPMAT(5X, 'MOVING  DOWN  USING  X( 

DC  123  KD=1,K 

B(KD)=B(KD)-X(LED,KD) 

WPITE(6,1001 )    KD,B(KD) 

CCNTINUE 

KEDS=8 

GC  TO  10 

PLOT  OUTPUT  ON  LINE  PRINTER 


' ,  13,  •  )  ' 


DC  825 

DO  827 

PICt  I ,  J 

CCNTINU 

CCNTINU 

DC  704 

PICQ,  J 

PIC(61t 

CCNTINU 

PIC(  1,1 

PIC(61, 

DC  705 

PIC(  1, J 

PIC(61, 

CCNTINU 

DO  714 

PIC(  I,  1 

PICU  ,3 

CONTINU 

PIC(31, 

PIC(31  , 

DO  715 

P  IC  ( 1 , 1 

PICt  I  ,3 

CONTINU 

DC  777 

JH=( YY2 

KH=(YY1 

IF(KH.G 

IFUH.G 

JF=38-J 

IF(PIC( 

IF(PIC( 

CCNTINU 

DO  888 

LC=(XX1 

LD=(XX2 

IFCLC.G 

IF(LD.G 

LD=38-L 

IF(PIC( 

IF(PIC( 

IF(PIC( 

IF(PIC( 

GC  TO  8 

PICCLC, 

CONTINU 

WRITE (6 

FCPMATt 

DC  895 


1=2 
J  =  2 
)  =  A 

E 

E 

J=l 

)=A 

J)  = 

E 

9)  = 

19) 

J  =  2 

)=A 

J)  = 

E 

1  =  1 

)  =  A 

7)  = 

E 

1)  = 

37) 

1  =  3 

)  =  A 

7)  = 

E 

JL  = 

(  JL 

(JL 

T.6 

T.3 

H 

KH, 

KH, 

E 

JN  = 

(JN 

( JN 

T.6 

T.2 

D 

LC, 

LC, 

LC, 

LC, 

88 

LD) 


,6C 
,36 
LPH(6) 


,18 

LPH(7) 

ALPH(7) 

ALPH(8) 

=ALPH(8) 

0,37 

LPH(7) 

ALPH(7) 

,30 

LPh(7) 
ALPH(7  ) 

ALPH(8) 

=ALPH(8) 

2,61 

LPH(7) 

ALPH(7) 

1,N 

)/SMAX)*17. 0+19.0 

)/TMAX)*29. 0+31.0 

C)  KH=60 

6)  JH=36 

JH).NE.ALPH(6))  PIC(KH,JH)=ALPF(2) 
JH)  .EQ.  ALPH(6) )  PIC(KH,JH)=ALPH( 1 ) 

1,M 

)/TMAX)*29. 0+21.0 
)/SMAX)*17. 0+19.0 
0)  LC=60 
6)  LD=36 


LD)  .E0.ALPH(6) ) 
LD)  .EQ.ALPHt  1) ) 
LD)  ,EQ.ALPH(2)) 
LD)  .EQ.ALPH( 3) ) 

=ALPH(3) 


GO  TO  890 

PIC(LC,LD)=ALPH(5) 
PIC(LC,LD)  =  ALPF(5  ) 
PIC(LC,LD)=ALPH(4) 


E 

,751) 
•1'  ) 
J=l,37 
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896 

895 


C 

C 
C 


68 
897 


1000 
10C1 
5O0C 
5001 
183 


WRITE(6,896)  (  P I C (  I  , J ) , 1= 1 , 61 ) 

FOPMAT(  20X,61A1) 
CONTINUE 
WPITE(6,751) 
GC  TO  55 

REGRESS  TO  THE  PREVIOUS  ITERATION 


WRITE(6,897) 


FCPMAT(5X 


EQ 
EQ 
EQ 
EQ 
EG 


IF(KEDS 

IF(KEDS 

IFCKEDS 

IFtKEDS 

IF(KEDS 

IF(KEDS.EQ 

IF(KEDS.EQ 

IFtKEDS. EC 

GO  TO  55 

F0PMAT(5X, 

FCPMAT(5X, 

FORMAT ( 5X, 

FORMAT (5X, 

STOP 

END 


REGRESSING  TO 

1)  GO  TO  60 

2)  GO  TO  59 

3)  GO  TO  62 

4)  GO  TO  61 

5)  GO  TO  64 

6)  GO  TO  63 

7)  GO  TO  66 

8)  GO  TO  65 


LAST  STEP' ) 


NEW  VALUE 
NEW  VALUE 
XXK  •  ,  13, 
YYK  '  ,  13, 


OF 
OF 
)  = 
)  = 


A( 
B( 
iF 

iF 


•f  13, 

•  ,13, 
10.4, 
10.4, 


•  )  =  • 

•  )  =  • 

»XX2( 
•YY2( ■ 


,F10 
,F10 


,13, 
,13, 


•  )  =  • 


,F10.4) 
,F10.4J 


PREPREPARED  SUBROUTINE  FOR  DATA  INPUT  AT  THE  AGT 


SLBROUTINE 

DIMENSION 

IE#IBLKS1 

NULL#-l 

IF3IBUF*l< 

1  IF3I6UF?1< 
ENC0DE316, 

10  FCPMAT^aNA 
CALL  TEXTO 
IF2IER .NE. 
RETURN 

50  CALL  TEXTO 
IFSIER.NE. 
CALL  TEXTO 
IFSIER.NE. 
RETURN 

100  CALL  TEXTR 
IF2IER.NE. 

110    IF2MOD3ITD 
CALL  TEXTI 
IFIIER.NE. 
RETURN 
END 


GINP2IDEV, ITDIR, IBLK, IEUF< 
IBUFSK,  ITDIRSK 


<1,50,  ICC 

•  NE.-KGO  TO  100 

10, IBUF< 

MEL  1ST  INPUT2< 

3IDEV,  IBUF,4,38,1 ,1,3,  IER< 

0<OUTPUT%10KIER,aGINPia 

%IDEV,NULL,1 ,38 ,1 ,1,3, IER< 

C<OUTPUT310KIER,3NULLla 
SIDEV, NULL, 1,40,4, 1, 3,  IER< 
0<OUTPUT%10KIER,3NULL23 

%  I DEV, NULL, 1,40,4, 1,3,  IER< 
0<OUTPUT%101<IER,3GINP2a 
IRS5IB<,8<.EQ.O<GO    TO     UC 
2IDEV,  IBUF,24,0,  IB,  IER< 
0<OUTPUT%10KIER,aGINP33 


$ 

C 
C 

c 

c 

-^ETA9300    SIfLOtGO 

$ 

$FNS 


PREPREPARED    PROGRAMS        IN    MACHINE    LANGUAGE    FOR 
THE    USE    OF    THE    LIGHTPEN    AND    THE    FUNCTION    SWITCHES 


IDEV 

ISV, 

IER 


PZE 
BRM 
PZE 
PZE 
PZE 
PZE 
LDA 
ERK 
BRU 
STA 
STA 


0 

9SETUPN 

3 

0 

0 

0 

*IDEV 

DVNOCK 

FNSDER 

FNSCLS2 

FICBFL 
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LDA 

#C10COOOO 

STA 

FNSSWO 

ECM 

032004 

FNSCL 

ERM 

D>EXEC 

PZE 

2 

PZE 

0 

PZE 

FNSSWO 

ECM 

022001 

ECM 

032002 

LDA 

*IDEV 

SKA 

IOBFL 

BRU 

$-1 

SKE 

#1 

BRU 

$£3 

LCA 

077766 

BRU 

$&2 

LCA 

077771 

STA 

*ISW 

LDA 

#C 

BRU 

$G2 

FNSDER 

LCA 

#1 

FNSER 

STA 

*IER 

BRR 

FNS 

FNSSkO 

PZE 

0 

PZE 

0 

PZE 

0 

END 

$ 

-iECF 

-.MET/S9 
* 

3CC  SI 

,LO,GO 

$GINFUT  PZE 

0 

BRM 

9SETUPN 

PZE 

3 

IDEV 

PZE 

0 

ITDIR 

PZE 

0 

IBLK 

PZE 

0 

LCP 

IDEV 

STD 

TEXT0G3 

LDA 

*IBLK 

STA 

BLOCK 

LDA 

INPADR 

ADD 

BUF 

STA 

IBUF 

LDA 

R>ONES 

STA 

*IBUF 

BRM 

TEXTO 

PZE 

0 

LCA 

INPADR 

ACD 

READ 

STA 

PATCH 

LDA 

BRM 

XNA 

*PATCH 

STA 

BRM 

LCA 

#101 

STA 

*IBUF 

INPADR 

EPM 

9INPUT 

STZ 

*IBUF 

BRM 

TEXTO 

PZE 

0 

LCA 

BRM 

XMA 

*PATCH 

STA 

BRM 

* 

* 

BRR 

GINPUT 

* 
TEXTC 

PZE 

0 

BRM 

GINP 

PZE 

4 

PZE 

0 

PZE 

0 

88 


PZE 

BLOCK 

IBUF 

PZE 

0 

NPO 

TEXTO 

* 

* 
BLCCK 

ERR 

TEXTO 

FZE 

0 

PATCH 

PZE 

0 

BR^ 

ERM 

TEXTO 

REAC 

DATA 

0563 

BUF 

DATA 
END 

0773 

$ 

-EOF 

-iLCAC 

XR 

-iDATA 
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APPENDIX  II .     PROGRAM  OPERATING  INSTRUCTIONS 

To  ready  the  program  for  a  set  of  data  a  total  of  eight  cards  must 
be  prepared  and  inserted  into  the  deck.     For  the  purposes  of  this  appendix, 
assume  that  there  are  GG  Y  points  and  HH  X  points.     Furthermore,  assume 
all  points  are  of  dimension  II.    Then  the  required  cards  are  as  listed 
below. 

Four  dimension  cards  are  placed  immediately  after  the  card  which 
reads  "DIMENSION  IGDIR(3)";  this  is  the  third  card  in  the  program.     These 
cards  should  read: 

DIMENSION  DAV(II),  TOT(II) 

DIMENSION  IMAG(6*HH+1)  ,IFRAM(5«GG+1)  ,A(II)  ,B(II) 

DIMENSION  X(HH ,  II) ,  Y(GG ,  II)  ,XX1(HH)  ,XX2  (HH) , 
YY1(GG),YY2(GG) 

DIMENSION  JMAG(6»HH+1),JFRAM(5*GG+1) 
Next  following  card  number  fifteen  of  the  program  which  reads 
"TT=0.  01" ,  the  cards  to  be  entered  are: 

M=HH 

N=GG 

K=II 
Finally  a  format  statement  must  be  inserted  into  the  program.    The  number 
of  the  statement  is  2187;  further  guidance  is  not  practical  since  the 
format  statement  is  used  in  connection  with  the  formatted  input  of  the 
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data.  As  an  example,  this  statement  might  read  "2187  FORMAT(4F10.  6) . 
The  statement  must  correspond  to  the  format  of  the  data  cards .  The  data 
cards  are  placed  behind  the  program,  and  the  deck  is  ready  to  run. 

Instructions  for  operating  the  XDS9300  computer  and  the  ADAGE  AGT 
console  are  contained  in  the  publication  entitled  "ELECTRICAL  ENGINEERING 
COMPUTER  LABORATORY"  by  R.  D.  Delaura  [2]   .     Familiarity  with  this 
document  is  a  must  before  operating  these  computers .    The  following 
instructions  are  a  guide  for  this  particular  program . 

1.  Energize  XDS9300 

a.  Depress  "RESET"  and  "POWER"  simultaneously,  and 
then  "SENSE  2"  on  the  9300  console. 

b.  Turn  the  teletypewriter  on. 

c.  Load  the  program  deck  into  the  card  reader  and  depress 
"POWER"  and  "START"  on  the  card  reader. 

d.  Depress  "READY"  on  the  line  printer. 

2.  Energize  the  AGT  Graphics  Console 

a.  Depress  "ON"  and  "RESET"  at  the  AGT  cabinet. 

b.  Turn  on  the  disk  drive. 

c.  When  the  ready  light  on  the  disk  cabinet  appears, 
energize  the  "THIS'IS  IT"  circuit  breaker  on  the  back 
of  the  AGT  cabinet . 

d.  On  the  AGT  cabinet  depress  "HALT",   "RESET",  "RUN", 
and  "PULSE  1". 
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e.  The  typewriter  at  the  graphics  console  should  type  the 
date  request:   "MO/DA/YR=".     Type  this  information  in. 
If  this  request  does  not  appear,  the  AGT  must  be 
"Bootstrapped".    Information  for  Bootstrapping  is 
available  in  the  laboratory. 

f.  Type  "RESET("MAD",pw) !  where  the  p  and  v  numbers 
indicate  the  location  of  MAD  on  the  disk.    PW  is 
normally  either  101  or  104.    When  the  carriage  of  the 
typewriter  returns  type  "GATED  !  "  . 

g.  Place  the  function  switch  overlay  shown  in  Figure  A-2-1 
over  the  function  switches  and  energize  the  lightpen. 
The  AGT  console  is  now  ready. 

3 .         Compile  and  execute  the  program 

a.  At  the  XDS9300  console  depress  "CLEAR"  and  "CLEAR 
FLAGS"  simultaneously,  and  then  "IDLE",   "RESET", 
"RUN",  and  "CARDS"      The  card  reader  should  start 
readying  cards,  and  the  line  printer  should  produce  a 
listing. 

b.  When  execution  commences,  the  teletype  should  print 
the  statement  "TYPE  IDEV=2*    AND  A  CARRIAGE  RETURN 
IF  USING  AGT  2".     The  input  light  will  then  come  on. 
Follow  the  instructions  if  using  AGT  2;  if  using  AGT  1, 
simply  type  an  asterisk  and  a  carriage  return.     The 
program  then  reads  the  data . 
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PROGRAM 

OVERLAY 

O 

O 

O 

PICK  Y 

O 

END  EDIT 

O 

RIGHT 

O 

LEFT 

O 

UP 

O 

DOWN 

O 

O 

o 

NEW  A  &  B 

O 

O 

PICK  X 

o 

PLOT 

O 

REGRESS 

o 

COMPLETE 
CHOICE 



Figure  A-2-1.     FUNCTION  SWITCH  OVERLAY 
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c.         At  the  AGT  console  the  word  "NAMELIST"  and  a 

blinking  cursor  will  appear  at  the  bottom  of  the  screen 
Type  in  the  initial  values  of  the  A  and  B  vectors .    Any 
accepted  vector  notation  is  acceptable;  the  following 
when  possible  (space  limitations  may  require  separate 
entry  of  component  values)  is  recommended: 

A=l ,  1 , 1 , 1  and  carriage  return 

B=l  ,1,1,1  and  carriage  return 

*  and  carriage  return 
If  any  mistakes  in  inputting  the  A  or  B  vectors,  the 
line  printer  will  output  an  error  message  when  the 
carriage  return  is  depressed.    Until  the  asterisk  is 
typed,  errors  can  be  corrected,  without  harm  to  the 
execution  of  the  program.    When  the  asterisk  and 
carriage  return  is  typed,  the  program  starts  computing 
projections  and  the  graphics  display  of  video  X's  and 
Y' s  appears . 
4.         Iterations  and  output 

a.         To  cause  an  iteration,  the  user  pushes  the  switch 

"PICK  X"  or  "PICK  Y"  as  appropriate.    Video  appears 
at  the  bottom  of  the  screen  indicating  that  the  program 
is  in  the  graphics  edit  mode.    The  user  depresses  the 
button  on  the  side  of  the  lightpen  and  places  the  tip 
of  the  pen  over  the  desired  point;  he  then  releases  the 
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button.     The  point  selection  is  indicated  when  a  portion 
of  the  picked  letter  vanishes.    When  the  desired  point 
is  so  indicated  the  operator  presses  the  function  switches 
"COMPLETE  CHOICE"  and  "END  EDIT".     Finally  he 
causes  the  iteration  to  occur  by  depressing  one  of  the 
switches:     "RIGHT",   "LEFT",   "UP",  or  "DOWN". 

b.  By  pressing  the  "REGRESS"  button,  the  user  can  negate 
the  effects  of  the  last  iteration. 

c.  If  the  operator  desires  to  initiate  a  new  series  of 
iterations,  he  depresses  the  "NEWA&B"  button.    The 
word  "NAMELIST"  and  the  flashing  cursor  on  the  screen 
and  the  program  is  effectively  at  step  3c. 

d.  Finally  by  depressing  the  switch  labelled  "PLOT",  the 
user  causes  the  line  printer  to  output  a  six  by  six  inch 
representation  of  the  situation  he  views  on  the  AGT 
screen. 

5.         Securing  the  computer 

a.         There  is  no  feature  in  the  program  which  causes 

execution  to  terminate.    Instead  when  the  operator  is 
finished  he  secures  the  computer.     First  in  securing 
the  AGT,  he  depresses  switch  "PULSE  1"  and  types 
"HOME!".     He  depresses  "HALT"  and  "RESET"  at  the 
AGT  cabinet  and  stops  the  disk  drive.     He  then  flips 
the  "THIS  IS  IT"  circuit  breaker  to  the  off  position. 
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b.  The  card  reader  and  the  teletype  are  turned  off. 

c.  The  9300  is  secured  by  depressing  "IDLE",  followed 
by  "RESET"  and  "POWER  ON"  (simultaneously). 
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APPENDIX  III .     DATA  LISTING 

This  appendix  provides  a  listing  of  the  data  sets  used  in  the 
illustrative  examples  presented  in  section  six  of  the  paper.    The  data 
is  formatted  so  that  for  each  case,  a  row  of  numbers  in  the  listing 
represents  a  single  tuple  for  that  problem.     Thus,  for  example,  the  data 
set  for  case  one  is  presented  in  three  columns,  while  that  for  case  six 
is  in  thirteen  columns.    As  an  aid  for  reference,  the  following  index  to 
this  appendix  is  provided: 

CASE  1  98 

CASE  2  102 

CASE  3  108 

CASE  4  120 

CASE  5  122 

CASE  6  12& 
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CAGE       1       DATA 


Yl 

Y2 

Y3 

It  4760 

1 

.0910 

1.6880 

1*5520 

1 

►  5980 

1.9390 

1.4770 

1 

.8920 

1.5140 

1.6330 

1 

.0480 

1.5200 

1.1660 

1 

.0040 

1.0620 

1.8390 

1 

.4130 

1.5910 

1.5960 

1 

.1130 

1.6650 

1.4980 

1 

.3760 

l.ono 

1.3660 

1 

.2080 

1.8670 

1.4940 

1. 

.5920 

1.7380 

1.5280 

1« 

.2360 

1.8400 

1.9900 

1. 

•  7610 

1.7300 

1.0630 

1« 

.8430 

1.7990 

1.2730 

1. 

.8040 

1.8820 

1.2510 

1« 

.7340 

1.9240 

1.4090 

1. 

.4700 

1.0110 

It  4660 

1. 

.5380 

1.4990 

1.7840 

1. 

.7510 

1.4130 

1.2930 

1. 

3150 

lt9900 

1.5360 

1« 

2590 

1.0840 

1.0230 

1. 

1780 

1.3750 

1.947C 

1. 

0900 

1.2740 

lt239C 

li 

1410 

1.5050 
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1.9490 
1.4020 
1.6640 
1.3350 
1.5510 
1.5730 
1.7510 


1 .5000 
1.2780 
1.1390 
1.3180 
1.0790 
1.1890 
1.1550 


1.7620 
1.6590 
1.6760 
1.4630 
1.4690 
1.9350 
1.9190 
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XI 

X2 

X3 

.1820 

.5940 

.7590 

.7510 

.4200 

.0270 

.3400       i 

.0070 

.0580 

.4540       , 

.0790 

.5710 

.3660       . 

.9420 

.7470 

.0100 

.3710 

.2370 

.0670 

.2880 

.3260 

.1290 

.1430 

.3530 

.3040 

.4830 

.2190 

.1950 

.7130 

.8210 

.0000       i 

.2570 

.6930 

.1160 

.1950 

.2770 

.7240 

,5050 

.3230 

.1800 

,8060 

.7180 

.7020 

,3840 

.3070 

.3450 

,3510 

.6360 

.5060 

,9670 

.2980 

.0290 

,1020 

.6560 

.7670 

,7870 

.5340 

•  1310 

.0920 

.1190 

.0380 

.8070 

.8260 

.0440 

.9880 

.2590 

.6430 

.5340 

.7950 

100 


.1000 

.9630 

.9820 

.2960 

.4110 

.1090 

.6920 

.8360 

.1050 

.3370 

.6720 

.0700 

.8820 

.4480 

.3050 

.3980 

.6590 

.5150 

.7680 

.3780 

.6270 
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CASE   2   DATA 


Yl 

1800 


5970 


1970 

5060 
5880 
0180 
5170 
7010 
1240 
8330 
6490 
0320 
9290 
0640 
0690 
C920 
5240 
1760 
5850 
6320 
1710 
6840 
2080 


Y2 

Y3 

1.7180 

1< 

►  1000 

1.6210 

1. 

.5260 

1.5970 

1. 

.5400 

1.5610 

1- 

.0350 

1.3280 

1« 

.6390 

1.0940 

1. 

.2890 

1.6270 

1. 

.3170 

1.4520 

1. 

.3420 

1.2390 

1. 

.7110 

1.1200 

2, 

.0000 

1.4200 

1. 

.9980 

1.3650 

1, 

.9130 

1.2240 

1, 

.8800 

1.9960 

1. 

.5410 

1.2910 

1. 

.9790 

1.2260 

1. 

.1740 

1.3370 

li 

.9190 

1.6900 

i< 

.4870 

1.6480 

1< 

.9620 

1.1750 

1. 

.6150 

1.6590 

1< 

.6430 

1.9320 

1< 

.4370 

1.2060 

1 

,5440 

Y4 


li 

.8920 

1. 

.3580 

1. 

.5900 

1* 

.1760 

li 

.6250 

1. 

.6460 

il 

.8120 

1. 

.4520 

1. 

.2150 

1. 

.0720 

1. 

.8610 

1. 

.2830 

li 

.7470 

1. 

,8370 

li 

,6380 

1. 

,8060 

li 

,1660 

li 

.6870 

li 

.2250 

li 

,5710 

li 

.2520 

1. 

.1620 

li 

.2670 

Y5 

1.6170 
1.5520 
1.7330 
1.8920 
1.7540 
1.2840 
1.1970 
1.8700 
1.5460 
1.3440 
1.4980 
1.4600 
1.4050 
1.0460 
1.1700 
1.2440 
1.5320 
1.4220 
1.3780 
1.7240 
1.8420 
1.3780 
1.4970 


Y6 

1.7030 
1.7460 
1.4840 
1.8990 
1.6840 
1.5200 
1.1370 
1.9170 
1.2090 
1.3970 
1.8100 
1.8250 
1.0800 
1.0390 
1.1950 
1.8020 
1.9480 
1.4130 
1.6270 
1.2160 
1.9580 
1.8270 
1.6530 
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1.9680 

1.3420 

1.5310 

1.4870      1 

.1530 

1.0740 

1.4900 

1.1330 

J. 5750 

1.9290      1 

.0840 

1.7150 

1.6920 

1.3280 

1.8990 

1.8300      1 

.1060 

1.2220 

1.5300 

1.9090 

1.5220 

1.2270      1 

.8050 

1.2390 

1.0030 

1.1580 

1.0800 

1.0960      1 

.7830 

1.4620 

1.3440 

1.0590 

1.7920 

1.2040      1 

.5150 

1.2900 

1.6490 

1.4  84  0 

1.6560 

1.2440      1 

.8740 

1.4870 

1.4610 

1.1810 

1.6360 

1.3970      1 

.9240 

1.1520 

1.7010 

1.3170 

1.7840 

1.4090      1 

.7470 

1.1710 

1-7450 

1.6550 

1.5510 

1.3890      1 

.3680 

1.8820 

1.3860 

1.0390 

1.6430 

1.2360      1 

.1580 

1.4720 

1.757C 

1.3710 

1.3980 

1.55C0      1 

.8610 

1.2130 

1.1280 

1.3960 

1.0450 

1.6470      1 

.2940 

1.4500 

1.2890 

1.6390 

1.7970 

1.6410      1 

.8250 

1.4830 

1.8480 

1.1710 

1.9780 

1.7740      1 

.5450 

1.6930 

1.5190 

1.0470 

1.6790 

1.4180      1 

.9150 

1.3880 

1.3120 

1.7410 

1.7030 

1.5280      1 

.0540 

1.7750 

1.7770 

1.0190 

1.4150 

1.2320      1 

.4420 

1.8390 

1.5040 

1.0720 

1.0800 

1.2360      1 

.5070 

1.2360 

1.4890 

1.2300 

1.2440 

1.3930      1 

.2870 

1.4500 

1.5990 

1.1760 

1.485C 

1.6160      1 

.4220 

1.6330 

1.9970 

1.1190 

1.4000 

1.3640      1 

.1440 

1.1940 

1.2390 

1.3550 

1.8640 

1.2330      1 

.0520 

1.6570 

1.8420 

1.1950 

1.9410 

1.1760      1 

.0520 

1.9500 

1.7240 

1.8400 

1.602C 

1.8330      1 

.9070 

1.7560 
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1.C070 

1.5280 

1 

,5400 

1 

,6480 

1< 

,5320 

1.3720 

1.9520 

1.9020 

li 

,7750 

1 

,2180 

1< 

,9300 

1.3440 

1.6680 

1.3920 

1. 

,6360 

.  1 

,5750 

1- 

,2920 

1.9370 

1.3950 

1.2630 

li 

,8380 

1« 

,0160 

1, 

,9050 

1.5380 

1.4920 

1.9890 

1. 

,5140 

1. 

,5000 

1« 

,6070 

1.8030 

1.6290 

1.3330 

1< 

•  0270 

1< 

,0400 

1< 

,7410 

1.2820 

1.9770 

1.6460 

1. 

,1450 

1« 

,9480 

1. 

,3720 

1.5850 

1 .^740 

1.9880 

1. 

,0880 

1« 

,5780 

1. 

,8820 

1.5570 

1 .9040 

1.3340 

1, 

•  1040 

1. 

,9740 

1. 

,0880 

1.0210 

1.1340 

1.7700 

1. 

,2590 

1. 

,1160  . 

1< 

,5510 

1.6540 

1.0290 

1 .2620 

li 

.1120 

li 

,5010 

1, 

,3440 

1.3500 

1.5820 

1.7010 

1. 

,5980 

1« 

,6560 

1, 

,0250 

1.6470 

1.5920 

1.4960 

1. 

,9630 

1« 

•  0200 

1« 

,8050 

1.7030 

1.2960 

1.1010 

1. 

,3690 

1. 

,0920 

1. 

,7890 

1.2970 

1.0590 

1.1460 

1. 

.5310 

1, 

,0580 

1. 

,9920 

1.6490 

1.1650 

1.1740 

1. 

,4850 

1. 

,7240 

1. 

,1400 

1.0050 

1.4300 

1.9070 

1. 

,0890 

1. 

,3600 

1. 

,1800 

1.0280 

1.7460 

1.9390 

li 

,3440 

1. 

,2960 

1. 

,6260 

1.26C0 

1.2030 

1.6010 

•  1, 

,5250 

li 

,6100 

1» 

,5750 

1.1510 

1.4310 

1.0500 

1< 

.2420 

1, 

,2700 

1. 

,1460 

1.8740 

1.2860 

1.1700 

1. 

9390 

li 

•  0110 

1. 

,4860 

1.7680 

1 »0390 

1.4970 

1. 

623C 

1« 

,9720 

It 

5320 

1.8760 
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XI 

X2 

X3 

X4 

X5 

X6 

.5230 

.2360 

.8400 

.9900 

.7610 

.7300 

.0630 

.3430 

.7990 

.2730 

.8040 

.8820 

•  2510 

.7340       « 

.9240 

.4090       . 

.4710 

.0110 

.4660 

♦  5380 

.4990 

.7840       . 

.7510 

.4130 

.2930 

.3150 

.9900 

.5360       . 

.2590 

.0840 

•  C230 

.1780 

.3750 

.9470       i 

.0900 

.2740 

.2390 

.1410       . 

.5050 

.9490       . 

.5000 

.7620 

.4020 

.2770 

.6590 

.6640       . 

.1390 

.6760 

.3350 

.3180       « 

.4630 

.5530 

.0790 

.4690 

.5730 

.1890 

.9350 

.7510       i 

.1550 

.9200 

.9710 

.9050       . 

.6570 

.1750 

.2750 

.6280 

.9490 

.4280       « 

.6740 

.4000       . 

.4000 

.1710 

.9670 

.2840       i 

.3050 

.4500       . 

.6880 

.0550 

.9230 

.4290       « 

.5430 

.4770       . 

.7090 

.6420 

.1120 

.1490       < 

.8840 

.5820      '« 

.2520 

.8250 

.6940 

.0800       . 

.4200 

.7150       . 

.0680 

.3000 

.1710 

.2650 

.3710 

.5280       « 

.4620 

.2820 

.5430 

.9090       < 

.6120 

.0640       < 

.5720 

.1500 

•  2330 

.9870       1 

.3460 

.3050       . 

.7480 

.8770 

.4390 

.3270 

.2310 

.7630       i 

.1980 

.8160 

.1820 

.5940 

.7590 

.7510       < 

.4200 

.0270 

♦  3340 

.0070       . 

.0580 

.4550 

.0790 

.5710 

.3660 

.9420       . 

.7470 

.0100       . 

.3710 

.2370 
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.0670 

.2880 

.3260 

.1290 

.1430 

.3530 

.3040 

.4830       i 

.2190 

.1950 

.7130 

.8210 

.0000      « 

.2570       , 

.6930       < 

.1160 

.1940 

.2780 

.7240 

.5050       . 

.3230       . 

.  18C0 

.8060 

.7180 

.7020 

.3840       < 

.3070 

.3450 

.3510 

.6360 

.5060 

.9670       . 

.2980 

.0290       « 

.1020 

.6560 

.7660 

.7870       < 

.5840 

.1310 

.0920 

.1190 

.0380 

.8070       < 

.8260 

.0440       i 

.9880 

.2590 

.6430 

.5340       . 

.7960 

.0980       . 

.9630 

.9820 

.2960       . 

4120 

.1090 

.6920 

.8360 

.1050 

.3370       i 

6720       i 

.0700 

.8820 

.4480 

.3050 

.3980       < 

.6590       » 

.5150       * 

.7680 

.3780 

.6270 

.4760 

.0910       . 

.6880 

.5520       . 

.8970 

.9890 

.4770 

.8920 

.5140 

.6330       » 

.0480 

.5200 

.1660 

.0040       . 

.0620 

.8390       < 

.4120 

.5910 

.5960       . 

.1120 

.6650 

.4980 

.3760 

.04J0 

.3660       « 

.2080 

.8670       « 

.4930       . 

.5920 

.7380 

.0600 

.1300       i 

.5420 

.2300       . 

.9750 

.7960 

.4330 

.3730       « 

.6580       < 

.4860       . 

.8540 

.4660 

.9550 

.5410       < 

.4120 

.8610       » 

.5130 

.5250 

.1120 

.3640       » 

.7990 

.8070       i 

.3800 

.0970 

.3250 

.3210       . 

.1200 

.0100 

.3380 

.9710 

.1340 

.6920 

.9610 

.1270 

.7700 

.1880 

.7670 

.7340       < 

.1320       < 

.4360 

.1450 

.4760 

•  23h0 

.6790 

.3150       < 

.6830       . 

.0800 

.9800 
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.3180 

.3060       « 

,5830 

,4550       i 

,1370 

6580 

.8900 

.4230 

,5380 

,1260 

1380 

,6260 

.6210 

.7120       « 

,4490       < 

,1480       t 

1080       i 

,2780 

.4210 

.1050 

,2450 

,2840       . 

,9340 

,7280 

.7180 

.1460       « 

i6220 

,7020       < 

,5230 

,5750 

.8260 

.6350       i 

,2660 

,8080       « 

,0580       « 

,5990 

.9400 

.1890       i 

,9330 

1780 

5400       i 

,4290 

.6360 

.9780       i 

,6290 

,3800       « 

,2920 

,8740 

.7800 

.2020       . 

,7030 

,0720 

,1500 

,8240 

•  9810 

.1380       . 

.6370 

,6800       1 

4160       « 

,8180 

.8840 

.2890       • 

,6940       « 

,3400       . 

8760       « 

,1090 

.7420 

.6210       « 

,9010       , 

1040       . 

,6380 

,5320 

.1650 

.0940       » 

,4250       , 

,4770       « 

,6980       . 

,2480 

.5580 

.5280       . 

,6880 

,2780       , 

,9530       < 

,2640 

.5700 

.1750       . 

,7170 

,9160       , 

,3200 

,2530 

.5110 

.0660       » 

.8640       . 

,4980       . 

,8720 

.1860 

.5590 

.0440       < 

,1750 

,7050       1 

,0160 

.5620 

.4030 

.9310 

,2990 

,6540       . 

.6000       < 

.7330 

.4630 

.4810 

.1160 

.0740       . 

.8700       , 

.2130 

.7390 

.0940 

.6510       « 

.9670       t 

.3490 

.1730 

.5040 

.6630 

.6380 

.7120 

,5130 

,8220 

.2120 

.0560 

,5920 

.6580       . 

.3910 

,5660 
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CASE.       3       DATA 


Yl 

Y2 

Y3 

Y4 

Y5 

Y6 

3. 

.1800 

3.  7 ISO       . 

.1000       t 

.8920 

.6170 

.7030 

2. 

5970 

3.6210       t 

.5260 

.3580-      « 

.5520 

.7460 

3. 

.1970 

3.5970       i 

.5400 

.5900       t 

.7330 

.4840 

3. 

5C6C 

3.5610       . 

.0350       t 

1760 

.8920 

.8990 

3. 

.5880 

3.3280       < 

.6390       , 

.6250 

.7540 

.6840 

3. 

C1S0 

3.0940       . 

.2890       t 

.6460       ■ 

.2840 

.5200 

3. 

5170 

3.6270       . 

.3170 

.8120 

.1970 

.1370 

3, 

•7010 

3.4520 

.3420       t 

.4520 

.8700 

.9170 

3, 

1240 

3.2390       i 

.7110 

2150 

.5460 

.2090 

3. 

8830 

3.1200      It 

.0000 

.0720       . 

.3440 

.3970 

3. 

6490 

3.4200       . 

.9980       , 

.8610 

.4980 

.8100 

3. 

•0320 

3.3650       « 

.9130 

.2830       . 

.4600 

.8250 

3. 

.9290 

3.2240       . 

.8800       t 

.7470 

.4050 

.0800 

3. 

.0640 

3.9960       . 

.5410       t 

.8370       ■ 

.0460 

.0390 

3. 

.0690 

3.2910 

.9790 

.6380       . 

.1700 

.1950 

3, 

.0920 

3.2260 

.1740 

.8060       . 

.2440 

.8020 
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.3600 

3, 

,1360       , 

,6270 

,2640       . 

.5870 

.7420 
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•  7300 

3.7550 

.8410 

3.1170 

.2940 

3.7700 

»662C 

3.4380 

.0^30 

3.6550 

.1760 

3.3640 

.4730 

3.5920 

.2480 

3.1010 

.5200 

3.8710 

.9160 

3.6580 

.5850 

3.6290 

.5820 

3.9750 

•  3700 

3.0850 

.4900 

3.6900 

.7370 

3.3540 

.4110 

3.3370 

.5760 

3.2860 

.8290 

3.2700 

.5800 

3.6950 

.3690 

3.2790 

.3820 

3.6450 

.1270 

3.1980 

.5720 

3.4330 

.2860 

3.5050 

.2860 

3.7920 

.3840 

.5200 

.5790 

.2400 

.8690 

.6100 

.9960 

.2700 

.2380 

.7960 

.7660 

.9990 

.0540 

.8750 

.8240 

.3610 

.6530 

.3430 

•  0160 

.7160 

.5420 

.0120 

.3760 

.3600 

.9530 

.7730 

.7360 

.9270 

.4220 

.4270 

.9370 

.8710 

.6050 

.4800 

.2220 

.0630 

.9100 

.7780 

.9610 

.2860 

.3620 

.9330 

.0310 

.1950 

.6110 

.1500 

.6590 

.5190 

.9250 

.9280 

.1540 
.0290 
.4730 
.7720 
.0340 
.7680 
.1360 
.1920 
.5650 
.0440 
.2680 
.2840 
.3100 
.2420 
.7840 
.1640 
.9120 
•  3430 
.6830 
.1270 
.2560 
.8560 
.0270 
.9240 
.6610 


.3580 

.7750 

.4320 

.3380 

.8290 

.0610 

.3120 

.3120 

•  0450 

.8730 

.7400 

.8930 

.1320 

.1100 

.9030 

.0440 

•  4140 

.5470 

.0340 

.4210 

.4850 

.8010 

.9090 

.7610 

.4420 
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.5830 

3.1080       . 

.2180 

.9280       < 

.2830 

.4910 

.4470 

3.9900       « 

.3380 

.2960       < 

.6000 

.2820 

.6280 

3.1680       . 

.0160 

.9830       . 

.0900 

.5360 

.0120 

3.5190 

.4070 

.6850       < 

.6030       . 

.3120 

.1260 

3.3290       « 

.7640 

.8650       « 

.7680 

.7830 

.0200 

3.2700       . 

.1700 

.1980       . 

.7710 

.8130 

.8310 

3.4950       < 

.5830 

.5500       • 

.2900       < 

.7670 

.7950 

3.1570       . 

.7200 

.5480       . 

.5000 

.0200 

.1740 

3.6320       . 

.6630 

.2370      . 

.8280 

.2300 

.4840 

3.7240       . 

.8260 

.9400      < 

.7660 

.4380 

.2480 

3.6520       « 

.2240 

.1180 

.3720 

.4540 

.3020 

3.5440       . 

.9520 

.9510 

.0800       « 

.5970 

.7180 

3.3850 

.0^00 

.9460       i 

.8180 

.6950 

.0330 

3.2690       < 

.7260 

.0820       i 

.9670 

.6790 

.8710 

3.457C       . 

.1220 

.4420       < 

.4670       . 

.0570 

.1810 

3.2470       < 

.1890 

.0560       . 

.9360 

.5260 

.3070 

3.0370 

.5560 

.7^00       < 

.3160 

.3080 
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CASE       4      DATA 


Yl 


Y2 


Y3 


Y4 


-  i 

68?0 

-  < 

.3590 

-.1320           -1. 

,8380 

V  , 

.7650 

1. 

,4370 

.6250  •             < 

,8080 

5220 

m   i 

,4680 

-.2980           -1< 

,8480 

-< 

.7580 

-. 

,5810 

1.7140                 . 

,3730 

1« 

,9940 

,0860 

-.1110 

,0550 

,7020 

•■  i 

,7940 

1.6710 

.2390 

3690 

-« 

,8250 

1.6220                . 

,5890 

»  1 

.8440 

-« 

•  2490 

.0820           -1. 

,7860 

*  ( 

.7320 

.0710 

-.4260           -1< 

,7740 

-» 

.7350 

,0330 

-.0700           -1. 

,8580 

1, 

.9130 

,3420 

.1820                . 

.2230 

-. 

,8180 

1. 

,4470 

-.9770 

.0940 

»  1 

,3800 

1« 

,4660 

-1.1350 

.6260 

••  4 

,5510 

1« 

,5800 

-.9870 

.4580 

1 

,9410 

«■ 

,3660 

.1630 

.0220 

,1130 

1 

,5680 

■sl.0630                 < 

.3450 

,0980 

1 

,5600 

.3780           -1« 

.1680 

1. 

,9060 

.2180 

-.4960                 t 
\ 

.1410 

1. 

•  0050 

*»  | 

.4050 

1.6320 

.1690 

1, 

,9770 

Wf 

•  290C 

.0450 

.0510 
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XI 

.1370 

X2 

►  6680 

X3 

-.5700 

X4 
.4500 

.2^+50 

.7640 

-.5690 

.0990 

•  1140 

.8740 

.2970       , 

.3620 

.4360 

.5580 

.A460       » 

.2870 

.1970 

.5910 

.4740 

.5380 

.0770 

•9320 

-.2750 

.1710 

.2170 

.6450 

-.6260       t 

►  1710 

.0970 

.5710 

-.7950       . 

.1120 

.2930 

.6520 

.2490       . 

.2060 

.2030 

.7060 

^.4710 

.3580 

.1990 

.8330 

.3580       . 

.1330 

.2670 

.7540 

.3010 

.2510 

.0280 

.9320 

-.3580      -« 

.0480 

.4860 

.7850 

.1990       i 

►  2100 

.2540 

.6090 

.2680 

.5980 

.1630 

.8880 

.2430 

.2220 

.4100 

.3490 

.3970       t 

.3460 

.2020 

.7230 

-.5960 

.1500 

.3550 

.5350 

-.6690       • 

.3320 

.2280       i 

.5960 

.4750 

.5370 
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CASE   5   DATA 


Yl 


Y2 


Y3 


Y4 


Y5 


Y6 


Y7 


Y8 


160 
155 
130 
161 
135 
165 
150 
148 
150 
120 
120 
190 
100 
150 

90 
120 
200 
120 
105 
210 

90 
110 
100 


.0 

5< 

.1 

u  < 

-0 

28. 

.0 

5< 

.2 

8« 

.0 

27. 

.0 

5< 

•  0 

C< 

-0 

26. 

.0 

5< 

.0 

.0 

27. 

•  0 

5. 

i0 

2. 

.0 

27. 

.0 

5. 

.0 

1. 

.0 

23. 

.0 

5. 

•  0 

9  < 

»0 

29. 

.0 

4. 

.8 

8< 

»0 

26. 

.0 

<H 

.5 

7. 

0 

31. 

.0 

4. 

.0 

6. 

.0 

33. 

.0 

5. 

►  1 

8« 

.0 

32. 

.0 

4. 

.2 

8. 

.0 

30. 

.0 

4. 

»4 

O  ( 
**   * 

.0 

35. 

.0 

4- 

•  0 

7, 

•  0 

29. 

.0 

4. 

.6 

8. 

.0 

30. 

.0 

4< 

.7 

7. 

.0 

35. 

4. 

.3 

.0 

30. 

.0 

4. 

.1 

8< 

.0 

28. 

.0 

5. 

•  0 

7. 

.0 

27. 

.0 

5 

.2 

8. 

»0 

26. 

.0 

4 

.0 

o. 

.0 

25. 

.0 

5 

.2 

1 

.0 

25. 

.0 

4 

.3 

.0 

25. 

0 

70.0 

450 

0 

85.0 

400 

0 

75.0 

560 

0 

70.0 

665 

0 

88.0 

570 

0 

95.0 

675 

0 

90.0 

580 

0 

85.0 

390 

0 

60.0 

435 

0 

55.0 

440 

0 

56.0 

650 

0 

55.0 

640 

0 

48.0 

430 

0 

65.0 

650 

0 

70.0 

655 

0 

67.0 

645 

0 

62.0 

660 

0 

63.0 

530 

0 

64,0 

435 

0 

67.0 

440 

0 

68.0 

430 

0 

60.0 

530 

0 

70.0 

440 

.0 

4. 

.10 

1. 

.0 

3. 

»12 

7. 

.0 

2« 

.96 

9. 

.0 

4. 

.04 

1. 

.0 

3. 

.05 

3. 

.0 

3« 

.72 

4. 

.0 

3- 

.83 

5« 

.0 

3< 

.75 

8. 

.0 

3. 

.05 

0. 

.0 

2. 

.78 

8. 

.  0 

2. 

.89 

3. 

.0 

3 

.46 

7. 

.0 

2. 

.98 

6. 

.0 

2. 

.84 

o. 

.0 

2. 

.96 

7. 

.0 

2 

.78 

C 

.0 

3 

.20 

3. 

.0 

3. 

.05 

9. 

.0 

3« 

.55 

7. 

.0 

3. 

.10 

0. 

.0 

3< 

.41 

1. 

.0 

2« 

.97 

5. 

.0 

3. 

.62 

7. 

122 


70 

100 

13C 

90 

80 

95 

70 

80 

95 

100 

95 

85 

70 

80 

85 

200 

260 

195 

195 

220 

220 

190 

285 

300 

225 


.0 

4 

►  5 

8 

.0 

23. 

.0 

4 

.8 

9« 

.0 

27. 

.0 

5. 

.2 

9. 

.0 

25. 

.0 

4. 

.5 

1. 

.0 

37. 

.0 

4. 

.6 

.0 

32. 

.0 

4 

.9 

0< 

•  0 

25. 

.0 

4« 

.4 

P, 

.0 

30. 

.0 

4. 

.2 

3. 

.0 

36. 

.0 

5. 

.1 

5 » 

.0 

31. 

•  0 

4. 

.6 

1. 

.0 

24. 

.0 

4. 

.S 

o. 

.0 

27. 

.0 

4. 

.7 

2. 

.0 

25. 

.0 

4. 

8 

1< 

.0 

26. 

.0 

5. 

.4 

0. 

-0 

21. 

.0 

5. 

.5 

3. 

.0 

33. 

.0 

3. 

.4 

C« 

0 

24. 

.0 

3. 

•  1 

8, 

.0 

21. 

.0 

3. 

»0 

9. 

.0 

20. 

.0 

3. 

.2 

9. 

.0 

19. 

.0 

3. 

.3 

0. 

.0 

24. 

.0 

3. 

.0 

8. 

.0 

25. 

.0 

3< 

.4 

.0 

26. 

#0 

3. 

.0 

l. 

.0 

19. 

.0 

3. 

.0 

o  . 
>  • 

.0 

20. 

.0 

3. 

.0 

0« 

.0 

22. 

64.3 

65.0 

70.0 

74.0 

78. 0 

82.0 

85.0 

69.0 

70.0 

76.0 

74.0 

73.0 

78.0 

80.0 

81,0 

98.0 

110.0 

105.0 

110.0 

95,0 

90. C 

96.0 

100.0 

102.0 

105,0 


450 
355 
380 
350 
450 
260 


262 


260 

265 

270 

355 

360 

365 

370 

355 

1210 

1220 

1130 

1010 

1205 

1210 

1070 

990 

1120 

985 


.0 

3« 

,40 

9. 

.0 

3. 

.10 

4. 

.0 

2. 

•  80 

9. 

,0 

2. 

.95 

1. 

.0 

2- 

.12 

2. 

,0 

3« 

.11 

3. 

,0 

3. 

.22 

4. 

,0 

3« 

.35 

6, 

♦  0 

3. 

.40 

9. 

,0 

3. 

.54 

0. 

.0 

3. 

.63 

1. 

,0 

3. 

.17 

9, 

.0 

3. 

.18 

0. 

.0 

2. 

.90 

2. 

.0 

3. 

.12 

5. 

,0 

4. 

.10 

o. 

,0 

3. 

.89 

6. 

,0 

2. 

.98 

4. 

,0 

3. 

.10 

6* 

,0 

3. 

.54 

9, 

.0 

3. 

.78 

4, 

,0 

3. 

.36 

o. 

,0 

4. 

.23 

1. 

,0 

3. 

.75 

0. 

.0 

3. 

.80 

5, 
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260 

280 

300 

310 

290- 

260 

290. 

160- 

2<f0 

195 

?90 

210 

180 

205 

215 

270 

290 

320 

210. 

210 

185 

200 

170 

140. 

110 


3 

3 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 
2 
3 
3 
3 
3 
3 
.3 
3 
3 
2 
3 
H 


0 
4 
0 

1 

0 
3 
1 
5 
1 
4 
5 
0 
9 
<+ 
1 
0 
2 
0 
0 
4 
2 
9 
0 
9 


0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 

c 

0 
0 
0 
0 
0 
0 
0 
0 
0 
0 


22.0 


20 
20 
19 
26 
22 
25 
20 
20 
21 
19 
22 
22 
23 
21 
20 
23 
19 
18 
21 
25 
26 
20 
20 
22 


97.0 

112.0 

108.0 

106.0 

94.0 

98.0 

100.0 

79.0 

88.0 

81.0 

94,0 

96.0 

97.0 

90.0 

100.0 

111.0 

102.0 

87.0 

112. C 

95.0 

96.0 

98.0 

95.0 

98.0 

130.0 


1090 

1200 

835 

1055 

1240 

1015 

1010 

1170 

990 

975 

860 

950 

990 

805 

700 

1170 

1350 

1160 

1010 

1190 

1055 

980 

1095 

990 

220 


.0 

2 

.94 

8. 

.0 

2 

.70 

1. 

.0 

3< 

.20 

5. 

•  V 

3< 

.45 

7. 

.0 

3. 

.50 

1. 

.0 

4. 

.60 

9. 

.0 

3. 

.56 

2. 

.0 

4. 

.26 

4. 

,0 

4. 

,40 

2. 

.0 

3. 

.58 

4. 

.0 

3. 

.40 

0. 

.0 

2< 

.70 

7. 

.0 

3. 

.83 

8. 

.0 

3. 

.25 

8. 

.0 

4. 

.05 

2. 

.0 

4. 

.60 

!• 

.0 

4. 

.42 

3. 

.0 

3. 

.64 

8. 

♦  0 

3. 

.81 

5. 

.0 

4. 

.00 

0. 

.0 

4< 

.18 

1. 

.0 

3. 

.46 

8. 

.0 

4< 

.15 

o. 

.0 

4< 

.25 

2. 

.0 

4. 

.51 

9. 

124 


100*0  5.6  8.0  19.0  128. C  216.0  2.60  4.0 

95.0  4.9  8.0  24.0  124. C  218.0  2.45  0.0 

65.0  6.2  9.0  30.0  134.0  200.0  2.61  2»0 

55.0  5.0  0.0  27.0  128.0  205.0  2.62  1.0 

70.0  5.3  7.0  28.0  118.0  204.0  2.35  2.C 

85.0  4.9  1.0  19.0  117.0  206.0  2.90  4.0 

115.0  5.0  0.0  21.0  122.0  198.0  3.25  7.0 

95.0  4.8  8.0  27. C  114.0  228.0  2.65  1.0 

120.0  6.1  9.0  24.0  120.0  244.0  2.82  8.0 
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CASE   6   DATA 


Yl 


Y2 


Y3 


Y4   Y5   Y6   Y7   Y8   Y9  Y10  Yll   Y12 


Y13 


1.   365. 

152 

It   517. 

126 

It   343. 

180 

It   432. 

110 

1.   376. 

115 

1.   410. 

110 

It   408. 

105 

1.   410. 

150 

It   315. 

118 

1.   546. 

150 

It   437. 

120 

It   506. 

106 

It   338. 

110 

It   446t 

120 

1.   426. 

135 

1 .   4 1 C  . 

140 

1.   483. 

120 

3.   377. 

114 

3.   410. 

130 

3.   414. 

120 

1.   369. 

140 

1.   335. 

132 

It   410. 

140 

90 

i   4. 

70 

.   It 

86« 

»   1. 

65. 

•   4t 

80. 

.   It 

7C< 

.   1. 

65. 

>   1. 

•   100. 

.   It 

76. 

»   1. 

70. 

.   4* 

7C 

>   4t 

t    72. 

4t 

70. 

•   2. 

70. 

»   It 

t    80. 

.   it 

85« 

.   1. 

70. 

.   1. 

70. 

»   2. 

85. 

.   2. 

80. 

>   4. 

so. 

»   4. 

32. 

.   3. 

80. 

.   3. 

2. 

4. 

3 

.   1« 

►   i. 

5. 

►   250. 

39t 

2. 

3. 

2 

►   2. 

»   2. 

2. 

►   326. 

38. 

1 1 

3. 

2 

►   2. 

>   3. 

2. 

.   27  8. 

42» 

1 » 

3. 

2< 

.   4. 

>   4. 

2. 

i   278. 

42. 

1  • 

3. 

2. 

.   1. 

►   3. 

2. 

.   278. 

41, 

1  . 

3. 

?« 

►   2. 

.'  3. 

2. 

.   255. 

43t 

1  • 

3. 

?• 

»   2. 

►   3. 

2. 

►   278. 

50. 

4. 

2» 

»   2. 

i   4. 

2. 

.   278. 

44, 

2. 

3. 

2. 

»   2. 

.   4. 

5. 

.   315t 

28. 

1. 

3. 

2. 

»   2. 

.   3. 

1. 

»   212. 

47. 

2. 

3. 

2. 

•   2. 

i   4t 

2. 

.   372. 

39, 

2. 

3. 

2. 

►   2. 

i   3. 

.   145. 

37. 

2. 

3. 

2. 

•   1. 

.   3. 

.   278. 

39. 

It 

2. 

2. 

•   1« 

.   4. 

.   278. 

50t 

It 

3. 

2. 

.   1. 

•   4t 

.   278. 

48t 

It 

3. 

2. 

.   2. 

.   4t 

»   250t 

35. 

2. 

3. 

2. 

•   2« 

.   4t 

»   278. 

48, 

2. 

2. 

2. 

•   2. 

•   4t 

5. 

.   350t 

38, 

It 

3, 

2. 

.   2« 

.   4. 

2. 

.   264. 

39, 

2. 

It 

2. 

.   i« 

■   3t 

5. 

.   246. 

49, 

It 

3t 

2. 

i   1. 

•   4* 

1. 

278. 

38, 

it 

3. 

i« 

.   1. 

•   4. 

1. 

.   278. 

38, 

it 

3. 

1. 

•   2. 

.   4. 

2. 

.   268. 

40, 

126 


1.   422 

.   152 

•    92 

»   1  .   2t 

It   517 

.   110 

80 

.   It   2. 

It   250 

►   130 

70 

.   It   2. 

1.   466 

i   114 

74 

It   366 

►   120 

90 

tt   271« 

►  no. 

70- 

I.   467. 

i   13C 

70. 

It   372. 

.   128. 

92. 

it   300' 

.   140« 

75« 

,   4,   2. 

It   348. 

•   110, 

80. 

It   288, 

•   130. 

80. 

It   477. 

.   115, 

so. 

It   406. 

.   140, 

•  no. 

.   2.   2. 

It   489. 
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